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ABSTRACT 

This  thesis  contains  the  analysis  of  data  of  message 
traffic  loads  on  the  Naval  Communications  Fleet  Multi- 
channel Broadcast  System  and  results  of  a  simulation  model 
of  the  system  under  various  channel  alignments.   Distri- 
butions of  interarrival  times,  message  lengtlis  and  requests 
for  screens  and  reruns  are  determined  from  the  data.   These 
distributions  are  used  in  a  simulation  model  of  the  Broad- 
cast System.   The  model  is  used  to  compare  the  average 
delays  caused  by  backlogs  of  messages  when  two  channels 
devoted  to  a  given  ship-type  are  a)  used  in  series  with  the 
second  channel  used  to  rebroadcast  messages  from  the  first 
channel  after  a  one  hour  delay,  and  b)  both  channels  are 
used  in  parallel  to  transmit  first-run  messages.   The 
results  of  the  simulation  show  that  backlogs  are  reduced 
considerably  by  running  the  two  channels  in  parallel  at 
all  times. 
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I.   THE  PROBLEM 

This  study  concerns  the  message  traffic  loads  on  the 
Naval  Communications  Fleet  Multi-Channel  Broadcast  System 
and  the  problem  of  determining  the  most  efficient  use  of 
the  8  traffic  channels  in  delivering  teletype  message 
traffic  to  the  ships  at  sea. 

The  Commander,  Naval  Communications  Command  message, 
R  021944Z  SEP  71,  officially  commissioned  this  study.   A 
copy  of  this  message  is  shown  in  Appendix  A  for  reference 
and  to  provide  an  example  of  a  Naval  Communications  message, 

Appendix  B  provides  background  information  on  the 
organization  of  the  Naval  Communication  System  and  a  des- 
cription of  the  Fleet  Multi-Channel  Broadcast  System. 

A.   THE  QUESTION  AS  POSED  FOR  SOLUTION 

The  original  question  as  it  was  posed  for  this  study 
was:   "When  a  significant  backlog  of  any  precedence  message 
exists  at  a  channel  pair  broadcast  position,  such  that  the 
primary  channel  delivering  first  run  traffic  is  overloaded, 
when  should  the  assigned  secondary  channel  be  activated  for 
first  run  traffic  and  when  that  backlog  ceases  to  be  signi- 
ficant, when  should  the  secondary  channel  be  deactivated 
from  first  run  traffic?". 

In  order  to  discuss  the  problem  further  we  must  first 
consider  how  the  primary  and  secondary  channel  pairs  are 
aligned  and  used  in  different  traffic  load  situations. 
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B.   TYPES  OF  CHANNEL  PAIR  OPERATION 

Table  XVII  in  Appendix  B  describes  the  f : iet  multi- 
channel broadcast  alignment  and  the  designated  use  for  each 
of  the  8  channels.   This  study  considers  only  6  of  the  8 
channels.   We  note  from  Table  XVII  that  these  6  channels 
form  3  pairs  of  channels  which  serve  different  types  of 
ships  copying  the  broadcast.   They  are:   channels  1  and  2 
for  the  Destroyer  force,  channels  3  and  4  for  the  Service, 
Amphibious  and  Mine  forces  and  channels  5  and  7  for  the 
major  warships.   Channels  2,  4  and  7  are  each  designated 
for  use  either  as  a  1  hour  delayed  broadcast  or  as  an  over- 
load channel  for  the  primary  channels  1,  3  and  5  respectively, 
This  allows  these  secondary  channels  to  be  employed  in  one 
of  two  different  ways  depending  on  the  traffic  load 
situation. 

1.   Situation  1 

When  the  traffic  load  for  a  channel  pair  is  such 
that  the  primary  channel  is  not  significantly  backlogged  the 
secondary  channel  is  used  for  the  1  hour  delayed  rebroadcast 
of  the  primary  channel's  traffic.   This  alignment  thereby 
provides  the  subscribers  of  the  channel  pair  two  chances  to 
copy  a  message  correctly  and  completely.   Thus  if  a  sub- 
scriber misses  a  message  or  a  portion  of  a  message  on  the 
primary  channel  the  subscriber  may  then  wait  for  one  hour 
and  attempt  to  copy  the  message  on  the  rebroadcast  channel. 
If  the  subscriber  should  miss  the  message  both  times  he  then 
must  use  the  procedure  described,  in  Appendix  B,  for  obtain- 
ing the  missed  message.   This  configuration  assumes  that  by 
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providing  this  redundancy  in  traffic  delivery  that  the 
number  of  service  requests  for  reruns  will  b   reduced  com- 
pared to  having  only  one  opportunity  to  copy  the  message. 
Further,  it  is  noted  that,  for  a  subscriber  to  originate  a 
screen  request  and  to  receive  a  reply  and/or  the  desired 
rerun,  the  experienced  normal  time  required  for  such  action 
is  in  excess  of  one  hour.   Hence  the  1  hour  delayed  rebroad- 
cast  appears  to  provide  better  service  provided  all  sub- 
scribers use  the  system  as  this  alignment  intends. 

On  the  other  hand  there  may  be  situations  in  which 
the  urgency  of  the  message  is  such  that  the  subscriber  may 
not  want  to  wait  and  gamble  on  the  1  hour  delay,  only  to 
possibly  miss  the  message  again,  or  perhaps  equipment  prob- 
lems may  prevent  a  subscriber  from  being  able  to  copy  the 
rebroadcast. 

2 .   Situation  2 

This  situation  provides  for  the  case  when  the 
channel  pair  traffic  loads  are  heavy,  and  significant  back- 
logs do  build  up  such  that  a  decision  has  to  be  made  by  the 
NAVCOMMSTA  to  activate  the  secondary  channel  to  transmit 
first  run  traffic  in  the  same  manner  as  the  primary  channel. 
In  this  case  both  channels  of  the  pair  are  used  to  clear 
the  backlog  of  traffic  that  has  accumulated  and  to  handle 
any  anticipated  increase  in  traffic  arrivals  which  may 
cause  any  increased  backlog.   At  the  same  time,  however, 
this  configuration  precludes  any  facility  for  rebroadcasting 
the  first  run  traffic  as  was  described  in  situation  1  above. 
It  is  also  assumed  that  this  configuration  will  also  prompt 
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more  screen  requests  to  be  initiated  by  the  subscribers 
since  it  intuitively  seems  reasonable  that  mere  messages 
may  be  missed  since  twice  as  much  first  run  traffic  is 
being  delivered  to  the  subscribers  of  the  channel  pair. 

In  cases  where  a  broadcast  channel  has  no  delayed 
rebroadcast  facility  use  is  made  of  what  is  called  "heading 
recaps."   Heading  recaps  are  a  periodic  rebroadcast  of  the 
headings  of  messages  previously  transmitted  on  that  partic- 
ular channel.   Although  the  broadcast  of  the  heading  recaps 
will  consume  some  of  the  valuable  times  needed  for  trans- 
mitting messages  which  are  backlogged  it  is  reasoned  that 
this  feature  will  have  an  effect  of  cutting  down  on  the 
number  of  screen  requests  which  were  prompted  because  the 
subscriber  was  trying  to  determine  only  if  a  missed  message 
was  of  concern. 

Figure  1  diagramatically  displays  the  above  des- 
cribed situations,  where  channels  1  and  2  were  the  chosen 
channel  pair  for  illustration. 

C.   DISCUSSION  OF  THE  PROBLEM 

Under  conditions  when  no  "significant"  backlog  of  any 
message  precedence  exists  the  channel  pair  operation  will 
operate  as  described  in  situation  1  above.   However  when  a 
"significant"  backlog  starts  to  build  up  it  may  then  be 
necessary  to  use  situation  2  in  order  to  make  a  more  timely 
and  efficient  delivery  of  the  traffic. 

The  problem  of  this  study  is  to  determine  what  is  meant 
by  a  "significant"  backlog  and  thus  determine  the  point  at 
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Figure  1.   Types  of  Channel  Pair  Operation. 
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which  the  channel  pair  is  switched  from  situation  1  to  2 ,  or 
vice  versa.   At  present  this  decision  is  mad   by  (i)  observ- 
ing how  the  system  has  behaved  in  the  past,  (ii)  by  observing 
the  system  in  its  present  state,  and  (iii)  by  using  experience 
to  predict  the  future  state. 

The  experienced  communicator  knows  that  the  traffic 
arrival  rate  at  the  broadcast  position  folio ws  a  diurnal 
fluctuation.   He  has  also  observed  when  the  channel  pair 
operation  was  changed  to  the  situation  2  alignment  that  the 
service  desk  after  a  time  delay  experienced  a  surge  of  screen 
requests  since  the  rebroadcast  channel  facility  was  discon- 
tinued.  This  surge  of  screen  requests  which  produce  message 
reruns  will  most  likely  contribute  further  to  the  backlog 
and  thus  the  backlog  will  get  worse  rather  than  better. 
Therefore  when  the  backlog  becomes  significant  the  decision 
maker  may  decide  to  hold  the  system  in  situation  1  and  hope 
that  things  get  better,  or  he  may  make  the  decision  to  go  to 
situation  2  and  hope  that  the  surge  of  screen  requests,  does 
not  contribute  to  the  backlog  in  a  significant  way. 

This  method  of  decision  making  leaves  something  to  be 
desired  and  this  study  was  conducted  in  order  to  investigate 
the  two  alternate  situations  and  to  determine  values  of  the 
decision  variable  for  when  switching  should  take  place. 

When  the  NAVCOMMSTA  activates  the  secondary  channel  for 
first  run  traffic,  rather  than  as  a  rebroadcast  channel, 
there  is  a  30  minute  delay  requirement  between  notification 
to  the  subscribers  and  the  time  the  channel  is  activated. 
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This  delay  time  used  by  NAVCOMMSTA  San  Francesco,  and  assumed 
common  among  other  NAVCOMMSTA * s ,  is  needed  for  subscribers 
to  prepare  necessary  equipment  to  copy  the  secondary  channel 
all  the  time.   Notification  of  the  intent  to  activate  the 
secondary  channel  is  sent  to  the  subscribers  on  the  primary 
channel. 
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II.   THE  APPROACH  TO  THE  PROBLEM  AND  ITS  SOLUTION 

A.   THE  DECISION  VARIABLE 

In  order  to  answer  the  original  question  it  is  first 
necessary  to  determine  when  a  significant  backlog  of  messages 
exists.   Let  us  say  that  a  backlog  is  considered  significant 
if  the  expected  or  average  time  a  message  of  a  certain 
precedence  spends  in  the  broadcast  queue  and  transmission 
system  exceeds  the  limits  of  the  speed  of  service  criteria 
shown  in  Table  XVI  in  Appendix  B. 

Actually  the  time  standard  referred  to  here  is  smaller 
than  the  speed  of  service  criteria  in  the  table /  because  the 
message  has  already  spent  some  time  in  transmission  before 
reaching  the  broadcast  position. 

During  the  analysis  of  screen  request  messages  in  this 
study  it  was  observed  that  these  messages  took  an  average 
of  2  hours  in  transmission  from  the  originator  to  the  ser- 
vice desk  in  the  NAVCOMMSTA.   If  we  observe  that  these 
screen  requests  are  almost  always  of  Immediate  or  Priority 
precedence  and  if  their  numbers  were  evenly  divided  then 
the  average  speed  of  service  criteria  would  be  approximately 
150  minutes.   Thus  80  percent  of  the  average  transmission 
time  was  used  in  transmitting  the  screen  request  messages 
from  the  ship  originators  to  the  NAVCOMMSTA.   If  for  example, 
these  messages  were  not  screen  requests  messages  destined 
for  the  service  desk,  but  rather  for  further  relay  on  the 
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broadcast,  we  can  see  that  on  the  average  these  messages 
could  spend  30  minutes  waiting  at  the  broadc  st  position 
for  transmission  before  exceeding  the  average  speed  of 
service  criteria.   We  allow  however  that  not  all  messages 
arriving  at  the  broadcast  position  are  from  ships  at  sea 
but  that  many  are  from  shore  establishments  and  hence  the 
messages  were  received  at  the  NAVCOMMSTA  via  the  AUTODIN 
network  which  has  a  much  faster  transmission  rate  than  the 
100  wpm  ship-shore  teletype  circuits.   Nevertheless  it  does 
not  appear  unreasonable  to  assume  that  an  arbitrary  but 
conservative  revised  speed  of  service  criteria  standard  to 
be  used  at  the  broadcast  could  be  50  percent  of  the  criteria 
shown  in  Table  XVI  in  Appendix  B. 

Reducing  the  times  in  Table  XVI  by  1/2  would  result  in 
times  shown  in  Table  I  below. 

TABLE  I 
BROADCAST  SPEED  OF  SERVICE  CRITERIA 


Precedence 


Average 
Speed  of  Service 


Limits 


Flash  (Z) 

Immediate  (O) 
Priority  (P) 
Routine  (R) 


As  fast  as 
humanly  possible 

15  minutes 


lh   hours 


3  hours 


Same 

15-30  minutes 

30  minutes- 3  hours 

1V-12  hours 
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The  decision  variable  of  the  problem  can  then  be  defined 
as  the  expected  time  a  message  of  a  certain  precedence 
spends  in  the  broadcast  queue  and  transmission  system,  i.e., 
average  time  in  queue  plus  the  average  time  in  transmission. 
In  addition,  a  measure  of  effectiveness  can  be  defined  to 
be  the  time  difference  between  the  expected  time  in  system 
and  the  broadcast  speed  of  service  criteria. 

It  would  then  be  logical  to  conclude  that  if  messages 
of  a  certain  precedence  spent  more  time  on  the  average  than 
is  allowed  by  the  broadcast  speed  of  service  criteria,  the 
broadcast  system  is  not  making  timely  and  efficient  delivery 
of  the  traffic.   Thus  if  this  condition  exists  for  any 
appreciable  period  it  should  be  necessary  to  use  both  the 
primary  and  secondary  channels  to  deliver  the  traffic. 

B.   MODELING  THE  PROBLEM 

The  type  of  model  which  best  fits  this  system  is  a 
queueing  model.   The  customers  are  the  messages  and  the 
servers  are  the  transmission  channels.   The  system  has  four 
priorities  or  precedences  of  messages  of  which  the  highest 
precedence  is  Flash.   The  Flash  message  is  allowed  to  pre- 
empt any  lower  precedence.   When  a  flash  message  preempts 
another  message  the  preempted  message  is  later  rerun  from 
its  beginning  after  the  flash  message  transmission  has  been 
completed.   A  flash  message  is  run  three  times  to  ensure 
receipt  by  the  addressee.   The  preemptive  method  described 
is  known  as  the  Preempt  Rerun  method  as  opposed  to  the 
Preempt  Resume  method  where  the  preempted  messages  will 
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continue  transmission  from  the  point  at  whicl-  it  was  inter- 
rupted. 

The  modeling  of  the  feedback  loop  of  this  system  is 
probably  the  most  challenging  because  of  its  various  time 
delays.   The  feedback  loop  referred  to  is  that  part  of  the 
system  in  which  screen  requests  are  made  by  the  subscribers 
and  from  which  the  message  reruns  are  generated.   This  feed- 
back loop  has  three  time  delays  associated  wi  th  it  plus  a 
message  grouping  effect  and  a  change  in  precedence  effect. 

The  first  time  delay  is  the  period  of  time  between  the 
actual  time  the  message  was  missed  and  the  time  at  which  a 
screen  request  was  sent  out  by  the  subscriber.   This  time 
delay  was  not  analyzed  due  to  the  complexities  of  measure- 
ment.  The  second  and  third  time  delays  are  the  time  in 
transmission  of  the  screen  request  to  the  NAVCOMMSTA  and 
the  time  required  to  process  the  screen  request  at  the 
NAVCOMMSTA  broadcast  service  desk. 

The  grouping  effect  is  caused  by  the  fact  that  screen 
requests  normally  require  a  number  of  messages  to  be 
screened  and  this  in  turn  results  in  the  rerun  messages 
being  delivered  to  the  broadcast  position  in  groups  of 
messages . 

The  change  in  precedence  effect  is  something  that  has 
not  been  previously  mentioned.   It  is  observed  that  the 
normal  precedence  assigned  to  a  screen  request  is  either 
Immediate  or  Priority  even  though  some  of  the  messages 
asked  for  may  be  of  Routine  precedence.   This  is  done  because 
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the  originator  of  the  screen' request  senses  a  greater 
urgency  in  obtaining  the  missed  messages  because  of  the 
above  described  time  delays.   This  upgrading  of  precedence 
then  may  require  the  NAVCOMMSTA  to  answer  the  screen  request 
at  a  higher  precedence  than  the  precedence  of  any  message 
referenced  in  the  request. 

C.  REQUIRED  DATA 

The  data  required  to  analyze  this  problem  falls  into 
three  areas.   First,  data  to  determine  the  type  of  message 
arrivals  at  the  broadcast  position.   Second,  the  message 
length  distribution  must  be  known.   Since  the  transmission 
rate  of  the  servers  is  constant  at  100  wpm,  the  message 
length  distributions  translate  into  the  service  time  dis- 
tributions.  Third,  data  concerning  the  screen  requests  is 
required.   It  is  essential  to  investigate  the  time  delay 
involved  from  the  request  for  a  rerun  until  it  is  received 
by  the  subscriber,  as  well  as  various  other  parameters  which 
will  be  more  fully  described  later. 

D.  SUMMARY  OF  RESULTS 

In  order  to  realistically  model  the  problem,  data  for 

interarrival  times  of  first  run  messages,  message  lengths, 

and  requests  for  sceens  and  reruns  was  analyzed.   The 

following  conclusions  were  made. 

i.   The  first  run  message  interarrival  times  were 
fitted  to  the  exponential  distribution  and  thus 
the  arrivals  follow  a  Poisson  Process.   (See 
Table  II.) 
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ii.   The  message  lengths  were  fitted  to  a  hyperexpo- 
nential  distribution  with  two  prr  portion 
parameters.   (See  Table  III.) 

iii.  The  interorigination  times  for  screen  requests 
fit  the  exponential  distribution  and  thus  were 
also  concluded  to  arrive  in  a  Poisson  Process. 
(See  Table  X.) 

iv.   The  transmission  delay  distribution  for  screen 
requests  was  fitted  to  a  right  shifted  expo- 
nential distribution.   (See  Table  X.) 

v.   The  delay  of  a  screen  request  at  the  service 
desk  was  concluded  to  follow  rig.it  shifted 
hyperexponential  distribution  for  high  load 
periods  and  a  right  shifted  exponential  distri- 
bution for  low  volume  traffic.   (See  Table  X.) 

vi.   The  number  of  rerun  messages  resulting  from  a 
screen  request  was  found  to  be  exponentially 
distributed.   (See  Table  X.) 

The  relationship  of  the  number  of  first  run  messages  and 
the  number  of  subscribers  with  the  number  of  screens  and 
reruns  was  analyzed  by  linear  regression.   The  results  indi- 
cate that  these  variables  are  independent  of  each  other. 
(See  Figures  12-19.) 

The  GPSS  simulation  model  of  a  broadcast  channel  pair 
was  used  to  simulate  the  FASW  channel  pair  in  various  con- 
figurations and  traffic  loads.   The  configurations  simulated 
and  compared  by  the  average  level  of  backlog  were:   i)  channel 
2  closed  to  first  run  traffic  and  used  as  a  rebroadcast  of 
channel  1  traffic,  ii)  channels  1  and  2  both  transmitting 
first  run  traffic  while  the  rebroadcast  facility  was  not 
used  and  iii)  channel  2  alternately  opened  and  closed  to 
transmitting  first  run  traffic  whenever  specified  backlogs 
were  exceeded.   The  result  of  these  simulations  was  that  back- 
logs were  reduced  considerably  by  running  the  two  channels 

in  parallel  transmitting  first  run  traffic  and  reruns. 
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III.   DATA  COLLECTION  AND  ANALY.  IS 

The  primary  source  of  data  used  in  this  thesis  was 
NAVCOMMSTA  San  Francisco  located  at  Stockton,  California. 
The  relative  close  location  of  Stockton  to  Monterey,  Cali- 
fornia allowed  frequent  visits  to  the  NAVCOMMSTA  both  to 
talk  to  personnel  at  the  station  and  to  collect  data. 

Although  the  fleet  atmosphere  of  the  Eastern  Pacific 
differs  from  other  parts  of  the  world,  it  is  felt  that  the 
type  of  data  collected  and  that  in  general  the  analysis 
results  are  typical  of  any  NAVCOMMSTA. 

All  data  which  was  analyzed  to  determine  a  goodness  of 
fit  to  a  theoretical  probability  distribution  was  analyzed 
by  using  the  Fortran  IV  computer  program  described  and  shown 
in  Appendix  C. 

A.   FIRST  RUN  MESSAGE  INTERARRIVAL  TIMES  ANALYSIS 

1.  Source  of  Data 

The  data  used  for  the  analysis  represented  4  days  of 
traffic  received  at  the  channel  1  broadcast  position  at 
NAVCOMMSTA  San  Francisco.   The  days  represented  were  5-8 
March  1971.   The  traffic  load  for  these  days  was  a  low  to 
medium  volume  of  142  to  401  messages  per  day. 

2 .  Analysis  Results 

In  order  to  show  that  the  messages  arrive  in  a 
Poisson  Process  it  is  necessary  to  show  that  the  interarrival 
times  are  distributed  according  to  the  exponential  prob- 
ability distribution. 
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The  data  was  analyzed  for  the  daily  aggregated 
arrivals  and  by  individual  precedences.   Tab- e  II  gives  the 
results  of  the  statistical  goodness  of  fit  tests.   We  see 
that  the  Kolmogorov-Smirnov  tests  indicated  rejection  of 
the  hypothesis  that  the  interarrival  times  are  distributed 
exponentially  while  the  Chi-square  tests  indicated  acceptance 
for  all  samples  except  one.   This  apparent  contradiction  is 
explained  by  noting  that  the  interarrival  times  were  measured 
to  the  nearest  integer  minute.   Hence  this  caused  the  data 
to  contain  groups  of  data  points  of  equal  value,  especially 
for  times  less  than  or  equal  to  3  minutes.   As  is  noted  in 
Appendix  C  this  grouping  effect  causes  the  Kolmogorov- 
Smirnov  (K-S)  statistic  to  inflate  and  hence  causes  rejection 

Figures  2  and  3  represent  the  graphic  results  for 
All  precedences  of  5  March,  the  data  for  which  both  statis- 
tical tests  indicated  rejection  of  the  hypothesis.   Figures 
4  and  5  represent  the  graphic  results  for  All  Precedences  of 
6  March  and  are  shown  for  comparison  to  the  5  March  graphs 
since  the  chi-square  test  indicated  acceptance  of  the 
hypothesis  for  the  6  March  All  Precedence  data. 

From  the  statistical  test  results  and  from  graphs  of 
the  data  it  was  concluded  that  the  interarrival  times  could 
be  reasonably  assumed  to  be  distributed  according  to  the 
exponential  distribution.   Since  this  data  represented  only 
channel  1  arrivals,  it  is  assumed  that  message  arrivals  at 
other  channel  pairs  also  arrive  in  a  Poisson  Process. 
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Figure  2.   Log  of  the  Tail  Distribution  for  Interarrival 
Times,  5  MAR.,  All  Precedences. 
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Figure  4.   Log  of  the  Tail  Distribution  for  Interarrival 
Times,  6  MAR.,  All  Precedences. 


29 


Interarrival  Times  (minutes) 
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6  MAR.,  All  Precedences. 
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B.   MESSAGE  LENGTH  ANALYSIS 
1.   Source  of  Data 

The  data  for  message  length  analysis  is  from  two 
sources.   The  first  source  contains  two  samples  of  messages 
and  was  taken  at  NAVCOMMSTA  San  Francisco.   These  samples 
represented  the  traffic  sent  on  5  October  19  71  on  channel  1 
(FASW)  and  channel  5  (FNSC) .   The  other  source  of  data  was 
the  Naval  Electronics  Laboratory  Center.   Th is  data  was 
obtained  from  the  World  Wide  Traffic  Demand   tudy  data  bank 
and  it  represents  all  messages  sent  on  a  fleet  broadcast 
in  various  parts  of  the  world  on  8  October  19  70.   This 
sample  therefore  represented  the  Fleet  Broadcast  System  as 
a  whole  and  could  be  assumed  not  to  be  biased  by  any  area 
effects . 

The  FASW  and  FNSC  samples  were  taken  by  measuring 
the  messages  in  teletype  tape  form.   This  method  of  sampling 
provides  for  a  high  degree  of  accuracy  in  determining  the 
length  of  a  message.   Teletype  tape  contains  100  characters 
per  foot  which  equals   20  words.   It  should  be  noted  that 
the  tape  includes  all  teletype  characters  which  are  processed 
by  the  transmitter  distributor  when  sending  a  message  such 
as  spaces,  blanks,  line  feeds,  shifts  to  upper  and  lower 
case  print,  etc. ,  characters  which  are  not  easily  identified 
on  a  printed  copy. 

The  message  lengths  of  the  World  Wide  sample  on  the 
other  hand  were  determined  by  measuring  page  copies  of  the 
messages  and  using  a  regression  equation  to  determine  the 
number  of  words  per  inch  of  page  copy. 
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2.   Analysis  of  Message  Lengths 

The  general  plan  for  this  analysis  w-  s  to  attempt 
to  fit  a  probability  distribution  to  the  message  length  data 
by  precedence  and  then  to  compare  that  fit  to  the  world  wide 
sample. 

The  first  fit  attempted  was  the  exponential  distri- 
bution.  The  results  of  this  showed  that  the  message  lengths 
were  not  distributed  exponentially. 

At  this  point  two  observations  of  tho  above  analysis 
were  made.   One  was  that  a  message  has  a  minimum  positive 
length,  since  each  message  has  to  have  a  heading  which  has 
a  date-time-group ,  an  originator  and  an  addressee  plus  a 
text.   It  was  estimated  that  such  a  minimum  would  be  approx- 
imately 20  words  long.   This  led  us  to  fit  the  data  with  a 
right  shifted  exponential.   However,  the  second  observation 
that  was  made  was  that  the  graphs  of  the  log  of  the  recipro- 
cal of  the  tail  distribution,  or  simply  called  the  log  of 
the  tail,  showed  that  the  data  in  every  case  had  a  decreasing 
failure  rate  at  least  for  message  lengths  up  to  approximately 
500  words.   One  particular  distribution  with  a  decreasing 
failure  rate  is  a  mixture  of  exponential  distributions  known 
as  the  hyperexponential  distribution. 

These  two  observations  led  us  to  fit  the  data  with 
a  right  shifted  hyperexponential  distribution  with  two 
parameters . 

The  remaining  question  was  whether  the  FASW  and  FNSC 
data  which  was  taken  from  one  day  only  was  representative  of 
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messages  sent  on  other  days  and  on  broadcasts  throughout 
the  world.   This  is  the  primary  reason  that  the  world  wide 
sample  was  obtained. 

It  was  found  that  the  world  wide  data  behaved  in 
the  same  manner  as  the  FASW  and  FNSC  data  and  that  the  right 
shifted  hyperexponential  made  a  very  good  fit.   It  was  noted 
however  that  one  of  the  major  differences  between  the  two 
sets  of  data  was  that  the  world  wide  data  had  very  few  data 
points  which  exceeded  1000  words  while  such  data  points  were 
very  much  in  presence  in  the  other  two  samples.   Since  the 
analyst  personally  took  the  FASW  and  FNSC  samples  it  was 
decided  that  these  large  data  points  were  not  freak  outliers , 
but  were  very  much  a  part  of  the  distribution.   On  the  other 
hand  it  is  noted  that  in  reality,  message  lengths  on  the 
broadcast  and  in  general  do  have  an  upper  bound  to  their 
length  which  is  estimated  to  be  approximately  2000  words. 
This  is  an  apparent  contradiction  between  the  empirical  and 
theoretical  distributions.   It  is  resolved  however  by  observ- 
ing that  the  probability  of  getting  such  a  long  message  is 
indeed  very  small. 

3.   Analysis  Results 

Table  III  shows  the  results  of  the  goodness  of  fit 
tests  for  message  lengths.   It  should  be  pointed  out  that 
for  all  goodness  of  fit  tests  that  not  only  were  the  statis- 
tical tests  used  in  making  a  judgment  on  the  fit  but  also 
interpretation  of  the  log  of  the  tail  and  the  Cumulative 
Distribution  Function  (CDF)  graphs.   As  is  noted  in  Appendix 
C/  the  Chi-square  statistic  for  large  sample  sizes  tends 
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to  inflate  and  careful  attention  must  be  made  to  the  selec- 
tion of  the  number  of  intervals  and  their  widths.   In  this 
analysis  these  drawbacks  of  the  chi-square  test  were  care- 
fully controlled.   When  the  results  of  the  K-S  test  and  the 
chi-square  test  conflict  we  choose  the  K-S  test  results 
unless  a  graphical  interpretation  could  explain  the  reason 
for  the  rejection  result. 

We  see  in  Table  III  that  the  K-S  test  results  indi- 
cate that  the  hypothesis  should  be  rejected  when  all 
precedences  are  aggregated  while  only  in  the  case  of  the 
FNSC  data  for  all  precedences  does  the  chi-square  test 
indicate  rejection.   Additionally  the  K-S  test  indicates 
rejection  for  Flash  and  Immediate  precedences  for  the  FNSC 
data.   It  is  noted  that  the  Flash  message  samples  are 
included  only  to  give  representation  of  the  parameters 
associated  with  this  precedence  since  the  respective  samples 
were  very  small  and  give  little  meaning  to  any  goodness  of 
fit  tests. 

Figures  6  through  11  show  the  log  of  the  tail  dis- 
tribution and  CDF  graphs  for  the  aggregated  precedence  for 
each  of  the  three  samples. 

In  the  cases  where  the  Table  statistical  results 
indicate  rejection  for  All  Precedences  of  each  of  the  three 
samples  we  see  that  the  empirical  curve  in  the  log  of  the 
tail  graph  deviated  from  the  theoretical  curve  in  the  message 
length  range  of  greater  than  600  words.   However  the  CDF 
graph  shows  a  reasonably  good  fit  to  the  theoretical 
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Figure  6.   Log  of  the  Tail  Distribution  for  FASW 
Message  Lengths,  All  Precedences. 
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Figure  7.   CDFs  for  FASW  Message  Lengths, 
All  Precedences. 
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Figure  8.    Log  of  the  Tail  Distribution  for  FNSC 
Message  Lengths,  All  Precedences. 
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Figure  9.   CDFs  for  FNSC  Message  Lengths, 
All  Precedences. 
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Figure  10.   Log  of  the  Tail  Distribution  for  World  Wide 
Message  Lengths,  All  Precedences. 
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Figure  11.   CDFs  for  World  Wide  Message  Lengths, 
All  Precedences. 
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distribution.   Since  the  hyperexponential  distribution  is  a 
mixture  of  two  exponential  populations  one  i  terpretation  is 
that  the  message  length  come  from  two  exponential  distribu- 
tions with  the  probability  or  proportion  parameters  p  and  1-p 
as  shown  in  Table  III.   It  seems  reasonable  that  all  message 
lengths  can  be  considered  to  be  distributed  according  to  the 
right  shifted  hyperexponential  distribution. 

C.   SCREEN  REQUEST  AND  FEEDBACK  LOOP  ANALYST  3 
1.   Reason  for  Analysis 

The  analysis  of  the  entire  screen  request  and  feed- 
back loop  of  this  queueing  system  is  a  key  part  to  the  solu- 
tion of  the  problem.   For  if  the  system  had  no  feedback  then 
the  solution  of  using  both  channels  of  a  pair  for  first  run 
traffic  is  obvious,  since  there  would  be  no  reason  for 
rebroadcasting  the  first  run  traffic. 

The  investigation  of  this  area  poses  several  new 
questions.   First  it  is  known  that  a  certain  percentage  of 
messages  are  missed  during  the  first  run  and  for  some  sub- 
scribers these  messages  are  successfully  copied  later  on  the 
rebroadcast  channel.   In  order  to  make  some  judgment  of  the 
effectiveness  of  the  rebroadcast  channel  we  can  ask,  "How 
many  of  the  subscribers  really  do  copy  the  secondary 
channel?"   Next  it  is  known  that  when  messages  are  missed, 
that  the  screen  requests  usually  contain  a  request  for 
screening  more  than  one  message.   This  implies  that  the 
subscribers  do  wait  for  a  period  of  time  before  they  ask 
for  screens  or  they  may  miss  several  messages  in  succession 
or  a  mixture  of  both.   This  time  delay  was  not  analyzed 
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since  it  would  be  very  difficult  to  collect  data  from  the 
subscribers . 

When  a  screen  request  is  originated  there  is  a 
transmission  delay  between  the  origination  time  and  the 
receipt  by  the  NAVCOMMSTA.   Also,  when  the  screen  requests 
are  received  by  the  NAVCOMMSTA  there  is  a  time  delay  for 
processing  the  request  which  is  related  to  the  number  of 
messages  that  have  to  be  screened  and  to  the  number  of  other 
screen  requests  waiting  at  the  service  desk. 

This  entire  process  results  in  the  feedback  of 
messages  to  the  broadcast.   The  distribution  of  rerun 
messages,  times  between  them,  and  their  relationships  to  the 
number  of  subscribers  and  first  run  messages  must  be 
analyzed. 

2 .  Source  of  Data 

The  data  for  all  of  the  screen  requests  and  feedback 
loop  analysis  was  obtained  from  NAVCOMMSTA  San  Francisco. 

3.  Number  of  Subscribers  that  Copy  Rebroadcast  Channel 
A  sample  of  14  days  of  data  was  available  for  this 

analysis.   During  the  period,  26  May-8  June  1971,  subscribers 
were  required  to  report  to  NAVCOMMSTA  San  Francisco  periodi- 
cally as  to  whether  they  used  the  rebroadcast  channel.   Addi- 
tionally when  a  subscriber  submitted  a  screen  request  he  was 
to  indicate  whether  he  was  using  the  rebroadcast  channel. 
Table  IV  indicates  the  results  in  average  numbers  for  each 
channel. 
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TABLE  IV 

NUMBER  OF  SUBSCRIBERS  USING 
THE  REBROADCAST  CHANNEL 


NR.  SUBS 
REQ  SCRN 
NR.  COPY   &  COPY 
NR.    NR.    NR.    NR.   NR.SUBS  SECONDARY  SECONDARY 
CHANNEL  1STRUN  SUBS  SCRNED  RERUN  REQSCRN   CHANNEL    CHANNEL 

FASW    293.1   14.6   26.2    6.5     3.9       3.4       0.8 

(23.4%)    (5.8%) 

FALD    243.4   19.5   24.8    5.4     2.4       2.6       0.4 

(13.5%)    (1.8%) 

FNSC    165.9   11.8   39.2    5.1     4.8       3.2       2.0 

(27.2%)   (16.9%) 


Although  this  is  a  small  sample  and  the  area  loca- 
tion for  the  sample  does  not  represent  the  system  as  a 
whole  it  does  however  show  that  the  secondary  channel  may 
not  be  used  to  its  fullest  benefit.   If  there  is  to  be  any 
conclusion  drawn  from  these  results  it  must  be  that  a 
surprisingly  small  percent  of  the  subscribers  indicated  that 
they  used  the  rebroadcast  channel.   An  explanation  for  this 
may  be  that  the  majority  of  the  subscribers  either  do  not 
desire  to  wait  for  an  hour  to  attempt  to  copy  the  missed 
message  and  instead  originate  a  screen  request  or  they 
obtain  the  missed  message  from  another  ship. 
4.   Analysis  of  the  Causes  of  Reruns 

For  this  analysis  two  sets  of  data  were  available. 
One  set  covered  the  period  of  26  May-8  June  1971  and  the 
other  the  period  of  25  August-17  September  1971.   The 
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information  in  the  sample  contained  daily  strtistics  on  the 
number  of  first  run  messages,  the  number  of  .  jbscribers 
copying  the  channel,  the  number  of  screen  requests  processed, 
the  number  of  messages  required  to  be  screened  and  the 
number  of  messages  which  had  to  be  rerun.   This  data  repre- 
sented four  different  channels,  FASW,  FALD,  FNSC  and  FRTT, 
the  single  channel  fleet  broadcast  which  is  not  part  of  the 
multi-channel  system. 

The  August-September  data  was  broken  into  two  parts. 
During  the  period  7-17  September  a  First  Fleet  Exercise 
named  ROPEVAL  3-71  was  conducted  and  during  this  period  the 
Operations  Evaluation  Group  (OEG)  conducted  an  experiment 
of  broadcast  channel  realignment  on  the  FMUL  broadcast  [1,2]. 
During  this  period  channel  4  was  removed  from  use  as  a 
secondary  channel  to  channel  3  and  was  converted  into  a 
dedicated  exercise  channel  designated  XRTT.   The  XRTT 
channel  carried  only  exercise  traffic  to  the  exercise 
participants  and  thereby  lightened  the  load  on  the  other 
channels  during  the  exercise.   This  experiment  had  no  direct 
relationship  to  this  study  except  that  the  data  collected 
during  the  pre-ROPEVAL  period  and  the  ROPEVAL  period  did 
provide  for  some  comparisons.   One  such  comparison  was  in 
the  channel  3-4  configuration.   Here  the  subscribers  had  a 
backup  channel  for  rebroadcast  for  one  period  while  during 
the  exercise  they  did  not  and  hence  this  would  hopefully 
provide  some  insight  into  the  effectiveness  of  the  rebroad- 
cast channel. 
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The  exercise  time  period  also  provided  for  moderately 
high  traffic  loads  and  a  larger  number  of  subscribers. 
Additionally  the  exercise  channel  XRTT  provided  an  oppor- 
tunity to  look  at  a  high  volume  of  traffic  being  transmitted 
under  simulated  wartime  conditions  while  this  channel  also 
had  no  backup  channel  for  rebroadcast  but  used  only  heading 
recaps  as  did  channel  3  during  this  period. 

Various  statistics  on  screened  and  rerun  messages, 
called  service  statistics,  are  shown  for  the  various 
channels  during  the  August-September  period  in  Tables  V-IX 
and  for  the  May-June  period  in  Tables  XVIII-XXI  in  Appendix 
D.   These  service  statistics  show  no  obvious  relationships. 

In  the  service  statistics  Tables  the  Service  Request 
Rate  is  abbreviated  SVC  REQ  RATE,  and  is  determined  by  the 
following  equation: 

SVC  REQ  RATE  =  ,Nr.  Messages  Screened 


(Nr.  First  Run) (Nr.  of  Subscribers) 


This  rate  is  used  by  some  analysts  of  communications  as  a 
measure  of  feedback  in  the  system  which  includes  the  rela- 
tionship of  first  run  messages  and  subscribers  [1] .   Due  to 
the  great  variation  of  this  rate  over  the  periods  analyzed 
it  was  not  clear  whether  there  is  in  fact  any  relationship 
among  these  three  measures.   Therefore  a  simple  linear 
regression  analysis  was  used  to  determine  if  any  such  rela- 
tionship existed  [3].   An  attempt  was  made  to  determine  if 
there  was  a  linear  relationship  between  the  following: 
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Number  of  First  Run  Messages  vs.  Number  of  Screens 
Number  of  First  Run  Messages  vs.  Number  of  Reruns 
Number  of  Subscribers  vs.  Number  of  Screens 
Number  of  Subscribers  vs.  Number  of  Reruns 
(Nr.  of  First  Run) (Nr.  of  Subscribers)  vs.  Nr.  of  Screens 
(Nr.  of  First  Run) (Nr.  of  Subscribers)  vs.  Nr.  of  Reruns 
These  relationships  were  chosen,  as  in  the  case  of 
the  SVC  REQ  RATE,  because  intuitively  one  would  suspect  that 
an  increase  in  the  number  of  first  run  messages  broadcast  in 
a  day  and/or  the  number  of  subscribers  copying  the  channel 
would  lead  to  an  increase  in  missed  messages  and  reruns. 

a.   Results  of  the  Statistical  and  Regression  Analysis 
An  important  item  in  Tables  V-IX,  and  those  in 
Appendix  D,  is  the  percentage  of  first  run  traffic  that 
result  in  reruns.   This  measure  indicates  the  additional 
work  that  is  placed  on  the  system  due  to  the  feedback.   We 
note  that  generally  this  percentage  is  between  2-4%  even 
though  the  number  of  first  run  messages  and  the  number  of 
subscribers  vary  from  channel  to  channel.   Two  exceptions 
are  8.2%  for  FNSC  (Aug-Sep)  and  18.2%  for  XRTT  during 
ROPEVAL  3-71.   Recall  that  ROPEVAL  was  an  exercise  simulating 
wartime  conditions  and  hence  this  increase  of  reruns  during 
this  period  is  not  surprising.   Probably  the  explanation  for 
this  is  that  the  subscribers  sensed  a  greater  urgency  in 
obtaining  missed  messages  plus  they  did  not  have  a  rebroad- 
cast  channel  associated  with  XRTT.   This  would  tend  to 
indicate  that  the  rebroadcast  channel  may  be  an  effective 
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feature  of  the  system.   On  the  other  hand  recall  that  during 
the  ROPEVAL  period  channel  3  was  operating  w  thout  a  rebroad- 
cast  channel  as  a  backup.   Looking  at  the  pre-ROPEVAL  and 
ROPEVAL  period  for  channel  3,  in  Table  VI  we  see  that  the 
percent  of  reruns  rose  1.3%  while  the  number  of  first  run 
messages  and  the  number  of  subscribers  were  also  greater 
during  the  ROPEVAL  period.   This  simple  comparison  appears 
to  indicate  that  the  presence  of  the  rebroadoast  channel  had 
little  effect  on  the  number  of  reruns.   An  important  appli- 
cation of  the  simulation  model  (described  later)  is  to 
resolve  these  apparently  conflicting  observations. 

Figures  12  through  19  show  representative  graphs 
of  the  regression  analysis.   All  graphs  represent  the  25 
August-17  September  period.   It  is  noted  that  the  regression 
lines  were  not  restricted  to  go  through  the  origin  even 
though  it  is  obvious  that  if  there  were  no  messages  trans- 
mitted there  would  be  no  feedback.   Therefore  the  reader 
must  be  cautioned  not  to  extrapolate  the  regression  line 
beyond  the  range  of  the  data  shown.   In  interpreting  these 
graphs  one  must  consider  the  grouping  or  spread  of  the  data 
points  and  the  coefficient  of  the  control  variable.   For  the 
FASW  graphs  in  Figures  12  and  13  we  see  that  the  regression 
line  is  nearly  horizontal,  i.e.,  it  has  a  slope  near  0. 
This  indicates  that  the  number  of  screens  and  reruns  are 
independent  of  the  number  of  first  run  messages  times  the 
number  of  subscribers. 
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The  channel  comparison  mentioned  above  is  shown 
in  Figures  14  and  15.   We  note  that  although  it  appears  that 
these  regression  lines  have  significant  slope  the  scale  of 
the  graph  causes  this  and  one  should  note  again  that  the 
coefficients  of  the  control  variable  are  nearly  zero  and 
hence  there  appears  to  be  no  significant  difference  in  the 
way  channel  3  was  operated. 

Figures  16  through  19  provide  results  for  the 
FNSC,  FRTT  and  XRTT  channels.   We  see  that  the  XRTT  regres- 
sion line  has  the  largest  slope  of  all  the  samples  shown, 
however  we  must  note  that  the  sample  size  is  nearly  one  half 
of  the  other  samples. 

If  any  conclusions  are  to  be  drawn  from  this 
analysis  it  would  be  that  the  number  of  reruns  are  indepen- 
dent of  the  number  of  first  run  messages  and  the  number  of 
subscribers.   However  further  analysis  in  this  area  appears 
necessary  in  order  to  fully  accept  this  conclusion. 

Regardless  of  what  conclusions  can  be  made, 
these  regression  equations  were  used  in  the  simulation  model 
(described  later)  to  mathematically  simulate  the  effect  of 
the  rebroadcast  channel  on  the  system. 

5.   Analysis  of  Time  Delays  in  the  Feedback  Loop 

Since  the  model  is  in  the  form  of  a  computer  simula- 
tion it  was  required  that  each  time  delay  of  the  feedback 
loop  be  analyzed  and  described  for  use  in  the  model. 

The  source  for  this  data  was  the  broadcast  service 
logs  from  NAVCOMMSTA  San  Francisco.   These  logs  contained 
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the  following  information:   the  Date-Time-Group  of  the  screen 
request  message,  the  time  the  screen  request  arrived  at  the 
service  desk,  the  time  the  screen  request  processing  was 
completed,  the  time  an  answer  to  the  screen  request  was  sent 
out  on  the  broadcast,  the  number  of  screens  each  request 
contained  and  the  number  of  reruns  which  resulted  from  each 
request.   Logs  from  two  different  months  were  used.   One 
represented  a  high  volume  traffic  month  which  was  September 
1971  and  the  other  a  low  volume  traffic  month  which  was 
December  1971. 

Four  time  distributions  were  analyzed.   These  were, 
the  interorigination  times  of  the  screen  requests,  the  inter- 
arrival  times  of  the  screen  requests  at  the  service  desk  in 
the  NAVCOMMSTA,  the  delay  in  transmission  times  for  the 
screen  request  messages  from  the  originator  to  the  NAVCOMMSTA 
and  the  delay  in  processing  the  screen  request  at  the  service 
desk.   Additionally  the  distribution  of  the  number  of  reruns 
per  screen  request  was  analyzed  for  use  in  the  model. 

a.   Interorigination  and  Interarrival  Distribution 
Analysis  Results 

For  each  of  these  distributions  the  times  between 
successive  origination  and  arrival  times  were  determined  and 
then  tested  for  goodness  of  fit  to  the  exponential  distribu- 
tion by  using  the  computer  analysis  program.   In  the  case  of 
the  origination  time,  the  date-time-group  of  the  screen 
request  was  considered  as  the  origination  time  of  the  request 
In  both  cases  the  hypothesis  that  these  times  were  dis- 
tributed exponentially  was  accepted.   The  results  are  shown 
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in  Table  X  and  Figure  20.   This  result  implies  that  screen 
requests  arrive  at  the  service  desk  accordir   to  the  Poisson 
Process  in  the  same  manner  as  first  run  traffic  at  the 
broadcast  position.   Only  the  interorigination  distribution 
is  used  in  the  model. 

b.  Transmission  Delay  Distribution  Analysis  Results 
The  delay  in  transmission  times  were  considered 

to  be  the  time  difference  between  the  origination  time  and 
the  time  of  arrival  at  the  broadcast  service  desk.   This 
data  was  tested  against  a  right  shifted  exponential  for 
goodness  of  fit.   The  right  shift  was  used  because  of  the 
observation  that  a  message  being  transmitted  from  a  ship  to 
the  NAVCOMMSTA  has  at  least  some  minimum  time  delay.   Table  X 
shows  the  results  for  both  the  regular  and  right  shifted 
exponential  fits.   The  hypothesis  that  the  data  fit  the 
right  shifted  exponential  was  accepted.   Figure  21  shows  the 
CDF  graph  for  the  September  data  sample.   It  is  interesting 
to  note  that  the  mean  transmission  delay  for  the  two  sets  of 
data  were  nearly  the  same.   The  sample  mean  for  the  December 
data  was  119.6  minutes  and  for  September  it  was  120.8 
minutes . 

c.  Service  Desk  Delay  Distribution  Analysis  Results 
The  service  desk  delay  times  were  considered  to 

be  the  time  difference  between  the  time  of  receipt  or 
arrival  of  the  screen  request  message  at  the  service  desk 
and  the  time  the  reply  to  the  screen  request  was  completed. 
It  should  be  noted  that  the  times  measured  here  included 
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both  the  time  the  screen  request  spent  in  qupue  at  the  desk 
and  the  time  for  processing  at  the  desk.   Thus  was  purposely 
done  in  order  that  this  distribution  could  be  used  in  the 
model  to  implicitly  simulate  the  delay  resulting  from  other 
screen  requests  for  other  broadcasts  and/or  channels  being 
handled  at  this  same  desk. 

This  data  was  tested  against  a  right  shifted 
distribution.   The  right  shift  was  again  used  by  the  obser- 
vation that  each  request  however  small  must  Lake  some  time 
to  be  processed.   Table  X  and  Figures  22  and  23  show  the 
results  of  these  tests.   In  the  case  of  the  December  data 
the  right  shifted  exponential  was  accepted  as  the  underlying 
distribution  while  for  the  September  data,  a  right  shifted 
hyperexponential  was  used  for  the  best  fit. 

d.   Reruns  per  Screen  Request  Distribution  Analysis 
Results 

The  data  for  this  analysis  represented  the 
number  of  reruns  per  screen  request  for  all  channels  pre- 
viously listed  for  the  August-September  period.   The  data 
was  fitted  with  a  regular  exponential  distribution  and  the 
hypothesis  was  accepted  that  the  reruns  per  screen  request 
data  is  distributed  exponentially.   The  results  of  the  test 
are  shown  in  Table  X. 
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Figure  20.   CDFs  for  Screen  Request  Interorigination 
Times  (September) . 
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Figure  21.   CDFs  for  Transmission  Delays  (September) . 
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Figure  22. 


Log  of  the  Tail  Distribution  for  Service 
Desk  System  Delays  (September) . 
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Figure  23.   CDFs  for  SErvice  Desk  System  Delays  (September) . 
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IV.       THE    MODEL 

A.  GENERAL  DESCRIPTION 

The  model  is  in  the  form  of  a  General  Purpose  Simulation 
System/360  (GPSS)  computer  program.   GPSS  is  a  simulation 
system  designed  for  simulation  of  queueing  systems  and  hence 
adapts  very  well  to  simulating  the  fleet  broadcast  system. 
References  4  and  5  fully  describe  GPSS. 

B.  MODEL  CAPABILITIES 

The  model  generates  transactions  which  are  the  messages 
in  the  system  and  moves  these  messages  through  the  model  for 
queueing  and  processing  in  a  very  similar  manner  as  the  real 
world  broadcast  system.   Since  each  of  the  channel  pairs  of 
the  broadcast  operate  independently  of  the  others,  the  model 
was  designed  to  simulate  only  one  channel  pair  of  the  fleet 
multi-channel  broadcast  system. 

The  model  has  essentially  3  basic  components  or  routines 
which  interact  with  each  other.   They  are  (i)  the  first  run 
traffic  flow  for  Immediate,  Priority  and  Routine  precedences, 
(ii)  the  Flash  message  routine  with  preemption,  and  (iii) 
the  screen  request  and  feedback  loop  routine.   For  the 
description  of  the  model  we  will  use  "channel  1"  to  refer  to 
the  primary  channel  and  "channel  2"  to  refer  to  the  secon- 
dary channel. 

The  model  is  essentially  very  simple  in  its  flow  charac- 
teristics and  has  much  flexibility  in  its  use.   This 
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flexibility  is  described  by  listing  the  following  features 
which  can  be  considered  as  measures  of  controlling  the  model 
The  user  may  change  or  specify  the  following  control  vari- 
ables, distributions  or  methods  listed. 

i.   The  daily  number  of  first  run  messages  to  be 
generated  for  the  simulation  period. 

ii.  The  percent  of  each  precedence  of  the  first  run 
messages  generated. 

iii.   The  individual  message  length  distributions  for 
each  precedence. 

iv.   The  mean  message  length  for  each  precedence. 

v.   The  diurnal  message  arrival  function. 

vi.   The  number  of  subscribers. 

vii.   The  backlog  levels  by  precedence  to  be  allowed 
before  opening  and  closing  channel  2  for  trans- 
mitting first  run  traffic. 

viii.   Allow  Flash  messages  to  preempt  other  messages 
using  either  the  Preempt  Rerun  method  or  the 
Preempt  Resume  method. 

ix.  The  mean  transmission  delay  for  screen  requests 
or  the  distribution  of  this  delay. 

x.   The  percent  of  screen  request  to  be  processed 

by  the  Broadcast  Control  Station  and  the  percent 
to  be  processed  at  other  NAVCOMMSTAs . 

xi.   The  service  desk  processing  delay  distribution. 

xii.   The  mean  service  desk  delay  for  the  BCS  and/or 
the  other  NAVCOMMSTAs. 

xiii.   The  mean  number  of  reruns  per  screen  request  or 
the  distribution. 

xiv.   The  mean  transmission  delay  for  reruns  to  reach 
the"  BCS  after  being  processed  by  the  other 
NAVCOMMSTAs . 

xv.  The  control  of  whether  rerun  messages  processed 
by  the  BCS  are  allowed  to  go  to  the  head  of  the 
line  of  first  run  messages  waiting  in  queue. 
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C.   DESCRIPTION  OF  THE  MESSAGE  FLOW  IN  THE  MODEL 

The  flow  charts  shown  in  Figures  2  4  through  2  7  provide  a 
simplified  picture  of  the  message  flow  through  the  model. 
First  consider  the  traffic  flow  for  Immediate,  Priority  and 
Routine  precedences.   These  messages  plus  Flash  are  generated 
according  to  an  exponential  distribution  with  a  mean  inter- 
arrival  rate  specified  and  a  specified  diurnal  function. 
After  a  message  is  generated  its  precedence  is  determined  by 
a  percentage  function  specified.   After  the  precedence  is 
assigned  the  model  routes  the  messages  by  precedence  to  be 
assigned  their  appropriate  message  lengths  from  functions 
specified  and  then  puts  the  messages  in  queue  by  precedence. 
At  this  point  in  the  model  the  Flash  messages  are  handled  in 
a  slightly  different  manner  than  the  other  messages.   This 
is  described  below.   The  precedences  or  priority  numbers 
assigned  by  the  model,  depending  on  the  precedence  assignment 
function  specified,  are  as  follows  (listed  from  lowest  to 
highest  with  designations  in  parentheses) 

1-Routine  First  Run  (R) 

2-Routine  Rerun  (RRR) 

3-Priority   First  Run    (P) 

4-Priority  Rerun  (PRR) 

5-Immediate  First  Run  (0) 

6-Immediate  Rerun  (ORR) 

7-Flash  (First  Run  and  Rerun)  (z) . 

The  model  chooses  the  highest  precedence  and  in  order  of 

arrival  within  precedences  when  advancing  transactions  through 

the  model.   Thus  messages  waiting  in  queue  are  chosen  by 

highest  precedence  and  on  a  first  in  first  out  basis  for 

routing  to  the  transmission  channels.   We  note  that  when 

priority  numbers  2,  4  and  6  are  assigned  to  reruns  processed 
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by  the  BCS,  these  messages  will  be  transmitted  before  their 
first  run  counterparts.   Hence  this  simulates  allowing  the 
reruns  to  be  put  at  the  head  of  the  line. 

When  a  message  is  routed  for  transmission  the  model 
checks  to  determine  the  number  of  backlog  by  precedence  and 
compares  it  to  a  specified  value  of  allowed  backlog  for  a 
given  precedence.   If  the  backlog  for  any  precedence  exceeds 
the  allowed  backlog  the  model  will  start  a  30  minute  timer 
(or  a  time  specified)  which  represents  the  notification  to 
the  subscribers  of  the  intent  to  activate  channel  2  for 
first  run  traffic.   At  the  end  of  the  30  minute  period 
channel  2  is  opened  for  first  run  traffic.   Whether  or  not 
the  timer  is  activated  the  message  proceeds  to  channel  1 
for  transmission.   Although  the  flow  chart  best  illustrates 
this  test  as  occurring  sequentially  between  the  time  the 
message  departs  the  queue  and  before  transmission,  this  is 
not  the  way  it  occurs  in  the  model.   Actually  all  the  above 
described  assignments  and  tests  are  made  at  the  instant  in 
simulation  time  when  the  message  is  generated.   That  is, 
when  the  message  is  generated  even  though  it  is  part  of  the 
queue,  the  message  proceeds  through  the  model  until  it  is 
refused  entry  to  some  GPSS  block  (program  statement) .   Thus 
the  messages  proceed  through  the  GPSS  model  to  the  BUFFER 
block  and  the  main  GATE  1  before  being  stopped.   For  statis- 
tical purposes  the  messages  do  not  depart  from  their  respec- 
tive queues  until  they  pass  through  the  DEPART  block  and 
this  occurs  only  when  a  channel  is  available  and  the  message 
is  allowed  to  SEIZE  a  channel  in  order  to  be  transmitted. 
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The  message  is  held  in  the  transmission  channel  (ADVANCE 
block)  for  a  period  of  time  equal  to  its  assigned  message 
length.   When  the  transmission  is  complete  the  message 
releases  the  channel  and  is  terminated  and  removed  from  the 
system.   The  main  GATE  1  closes  when  a  message  SEIZES  a 
channel  and  opens  when  the  message  RELEASES  a  channel. 

If  channel  2  has  been  activated  for  first  run  traffic, 
then  the  model  processes  the  messages  through  channel  2  in 
the  same  manner  as  in  the  primary  channel.   Messages  waiting 
for  transmission  are  selected  in  the  precedence  order 
described  above  and  each  time  a  channel  is  finished  with  a 
transmission,  or  is  available,  the  model  will  move  the 
available  message  to  the  appropriate  channel.   When  a  mes- 
sage is  processed  by  channel  2 ,  the  backlog  of  the  Priority 
and  Routine  precedences  are  checked  against  a  specified  value 
to  determine  if  channel  2  should  be  closed.   If  channel  2 
can  be  closed  this  will  activate  a  timer  which  runs  for  the 
duration  of  the  message  length  of  the  message  in  channel  2 
and  then  closes  channel  2  to  any  further  first  run  traffic. 

It  is  noted  with  emphasis  that  the  secondary  channel  in 
the  model  does  not  process  messages  when  it  is  considered  as 
being  used  for  rebroadcast  of  the  primary  channel  traffic. 
In  the  model  this  feature  is  handled  through  a  mathematical 
relationship  which  will  be  more  fully  described  below. 

The  next  part  of  the  model  to  be  described  is  the  Flash 
message  routine.  This  is  described  separately  to  point  out 
that  although  the  Flash  messages  are  generated  as  a  part  of 
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Figure  24.   First  Run  Traffic  Flow  in  the  Model 
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the  overall  first  run  traffic  described  above,  these  messages 
follow  a  different  route  through  the  model  because  they  are 
allowed  to  preempt  other  lower  precedence  messages.   When  a 
Flash  message  arrives  in  the  system  it  immediately  preempts 
any  message  in  channel  1  and  begins  transmission.   Just 
before  transmission  and  in  the  simulation  time  instant  of 
generation  the  Flash  message  is  assigned  its  message  length 
along  with  a   25  word  header  and  then  the  entire  length  is 
multiplied  by  three.   This  part  of  the  simulation  represents 
two  things.   First  before  a  Flash  message  starts  transmission 
on  the  broadcast  channel  a  heading  is  sent  which  reads, 

"ZUJ  FOR  FLASH • ZUJ  FOR  FLASH"  etc. ,  and  means  Stand 

By  for  a  Flash  Message.   This  header  also  causes  attention 
bells  to  ring  on  the  subscriber's  teletype  thereby  drawing 
attention  to  the  fact  that  a  Flash  message  is  about  to  be 
transmitted.   The  second  item  is  that  a  Flash  message  is 
always  transmitted  three  times  in  succession  to  ensure  re- 
ceipt by  the  subscribers. 

The  Flash  routine  also  provides  for  the  situation  when 
more  than  one  Flash  message  may  arrive.   In  this  situation 
if  channel  2  is  being  used  for  first  run  traffic  and  chan- 
nel 1  is  already  processing  a  Flash  message  then  the  second 
Flash  will  preempt  any  message  in  channel  2. 

When  a  message  is  preempted  by  a  Flash  message,  the 
preempted  message  is  handled  by  the  model  in  such  a  manner 
that  when  the  Flash  finishes  transmission,  the  preempted 
message  will  be  run  again  from  its  beginning,  if  the  model 
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is  being  operated  in  the  Preempt  Rerun  mode.   If  channel  2 
is  open  for  first  run  traffic  and  a  message  is  preempted  in 
channel  1  and  if  channel  2  becomes  available  before  the 
Flash  has  finished  transmission  in  channel  1  the  preempted 
message  will  be  routed  and  transmitted  on  channel  2. 

The  third  part  of  the  model  is  the  Screen  Request  and 
Feedback  Loop  of  the  broadcast  system.   Screen  requests  are 
generated  according  to  an  exponential  interorigination  dis- 
tribution with  specified  mean  value.   When  the  screen 
requests  are  generated  they  are  assigned  a  transmission  time 
from  a  right  shifted  exponential  distribution  with  a  given 
mean  time.   This  holds  the  request  in  an  ADVANCE  block  for 
this  period  of  time  before  allowing  it  to  proceed  through 
the  model.   It  is  noted  that  more  than  one  request  can  be 
in  this  ADVANCE  block  at  one  time  thus  simulating  any  number 
of  screen  requests  in  transmission  to  the  NAVCOMMSTA.   The 
model  then  randomly  determines  according  to  a  percent  speci- 
fied whether  the  screen  request  will  be  handled  by  the  BCS 
or  by  another  NAVCOMMSTA.   Next  the  request  is  assigned  a 
service  desk  system  processing  time  and  is  held  in  an  ADVANCE 
block  to  simulate  the  time  a  screen  request  spends  in  the 
service  desk  queue  and  on  the  service  desk  being  processed. 
This  time  delay  distribution  with  its  specified  mean  repre- 
sents the  competition  a  screen  request  for  the  channel  pair 
being  simulated  has  with  other  screen  requests  for  other 
broadcast  channel  pairs  or  other  broadcasts.   When  the 
request  has  finished  processing  and  is  ready  to  proceed, 
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it  is  then  assigned  the  number  of  reruns  that  have  resulted 
from  this  screen  request  from  an  exponential  distribution. 
If  the  number  of  reruns  is  equal  to  zero  the  screen  request 
transaction  is  terminated  and  removed  from  the  system.   If 
the  number  of  reruns  is  greater  than  or  equal  to  one,  then 
each  rerun  becomes  a  message  itself  and  is  assigned  a  prece- 
dence and  an  appropriate  message  length  and  sent  to  its 
appropriate  queue  to  await  transmission.   The  option  of 
allowing  reruns  to  go  to  the  head  of  the  line  was  described 
above.   If  a  screen  request  was  processed  by  a  NAVCOMMSTA 
other  than  the  BCS,  those  reruns  are  assigned  a  transmission 
time  from  the  transmission  delay  distribution  with  a  speci- 
fied mean  and  held  for  that  period  in  an  ADVANCE  block  before 
joining  their  respective  queues.   This  time  simulates  the 
time  it  takes  for  the  reruns  to  be  transmitted  via  AUTODIN 
to  the  BCS.   It  is  noted  that  in  this  case  the  reruns  cannot 
be  allowed  to  go  to  the  head  of  the  line  since  these  rerun 
messages  arrive  at  the  broadcast  position  in  the  same  manner 
as  first  run  messages  and  therefore  are  handled  in  the  same 
manner. 

D.   THE  MODEL  AND  THE  REBROADCAST  CHANNEL  SIMULATION 

As  was  mentioned  above,  this  model  does  not  rebroadcast 
first  run  messages  as  though  they  were  being  sent  on  a  1 
hour  delay  rebroadcast  channel.   The  secondary  channel  in 
the  model  is  used  solely  for  processing  first  run  messages 
and  when  it  is  not  doing  so  it  is  not  used.   The  effect  of 
the  use  or  non-availability  of  the  rebroadcast  channel  to 
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the  subscribers  is  simulated  through  the  control  of  the  mean 
interorigination  times  of  the  screens  requests.   This  is 
done  by  using  a  representative  regression  equation  of  (Nr. 
First  Run) (Nr.  Subscribers)  vs.  Nr.  Screens  for  the  channel 
pair  being  simulated.   Thus  by  specifying  the  number  of 
first  run  messages  and  the  number  of  subscribers  for  the 
simulation  period/  the  model  determines  the  number  of 
screens  that  will  be  handled.   Then  by  further  specifying 
the  average  number  of  screens  per  request  the  model  deter- 
mines the  number  of  screen  requests  to  be  generated  and  with 
this  information  the  model  calculates  the  mean  interorigina- 
tion time  by  dividing  the  number  of  requests  into  the 
simulation  time  period.   In  the  real  world  system  it  is 
assumed  that  when  the  rebroadcast  channel  is  being  used, 
that  the  number  of  screens  requests  is  less  than  when  the 
secondary  channel  is  not  being  used  for  rebroadcast  of 
first  run  traffic.   Then  if  this  assumption  is  in  fact  true 
the  analyst  need  only  decrease  the  mean  interorigination 
time  for  the  screen  requests  when  the  secondary  channel  is 
open  for  first  run  traffic.   This  therefore  simulates  a 
higher  number  of  screen  requests  being  generated  and  handled 
at  the  service  desk  and  thus  more  reruns  are  generated. 
This  is  done  in  the  model  by  one  of  two  ways.   The  regres- 
sion equation  used  can  be  changed  when  channel  2  is  open  or 
simply  the  number  of  screens  per  request  can  be  changed  in 
order  to  effect  the  increase  of  screen  requests.   Note  that 
this  change  in  the  mean  interorigination  time  is  made  at 
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the  time  that  the  secondary  channel  is  opened  or  closed  for 
first  run  traffic.   Thus  we  can  observe  how  the  increase  in 
the  number  of  reruns  affect  the  system  and  since  the  time 
delays  are  built  in  the  feedback  loop  we  can  also  observe 
any  surge  in  reruns  which  may  occur  after  the  secondary 
channel  was  opened  for  first  run  traffic. 

E.   VERIFICATION  OF  THE  MODEL 

As  the  model  was  built,  each  routine  added  was  run 
separately  to  verify  that  the  logic  design  was  correct  and 
the  transaction  flow  was  being  performed  correctly.   The 
means  and  variances  of  interarrival  times  and  message 
lengths,  etc.,  generated  by  the  model  using  random  number 
generators  and  the  specified  functions,  were  checked  for 
accuracy.   Appropriate  modifiers  (explained  in  model  docu- 
mentation) were  chosen  to  ensure  that  the  model  generated 
the  function  values  with  a  reasonable  accuracy  for  the  mean 
values  desired. 

One  important  verification  was  made  by  running  the  model 
for  the  M/M/l  and  M/M/2  queueing  systems.   That  is,  the 
model  was  run  without  the  feedback  loop  and  the  Flash  pre- 
emption routine  and  with  exponential  (M)  interarrival  and 
service  times  for  one  server  (only  channel  1  open)  and  two 
servers  (both  channels  open)  [6,7].   The  results  of  these 
model  runs  were  compared  to  analytical  calculations  for  the 
average  server  utilization  and  the  average  queue  length  and 
waiting  times.   These  comparisons  showed  that  the  model  was 
very  accurate. 
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The  model  was  also  verified  in  its  final  form  by  using 
data  obtained  from  OEG  which  was  gathered  during  ROPEVAL 
3-71   [1] .   This  data  represented  the  actual  hourly  number 
of  arrivals  and  backlogs  observed  during  the  period  of 
13-16  September  on  channel  1  (FASW) .   For  the  comparison 
runs  with  the  model,  parameters  of  the  model  were  set  at 
values  that  were  observed  during  this  time  period.   The 
parameters  which  were  not  varied  for  the  comparison  runs 
and  for  later  described  runs  are  shown  in  Table  XI  below. 

The  observed  traffic  loads,  number  of  subscribers  number 
of  screen  requests  and  number  of  reruns  are  shown  in  Table  XII 

TABLE  XI 
MODEL  PARAMETERS  HELD  CONSTANT 
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TABLE  XII 
FASW  13-16  SEP  OBSERVED  TRAFFIC  LOADS 


Nr. 

Nr. 

Nr. 

Nr. 

Screen 

Msgs 

Day 

First  Run 

Sub 

scribers 

Requests 

Rerun 

13  SEP 

422 

38 

6 

6 

14  SEP 

668 

40 

8 

4 

15  SEP 

598 

40 

13 

10 

16  SEP 

480 

33 

6 

36 

Totals 

2168 

33 

56 

To  simulate  the  daily  and  hourly  fluctuations  of  traffic 
arrivals  during  the  period  13-16  September  the  model  used 
the  number  of  first  run  messages  shown  in  Table  XII  above 
and  the  "diurnal"  function  shown  in  Figure  2  8  below.   The 
specified  dirunal  function  allows  the  model  to  determine 
the  proportion  of  a  day's  first  run  traffic  that  will  arrive 
or  be  generated  each  hour  of  the  simulation  run.   Another 
function  in  the  model  specifies  the  number  of  subscribers 
copying  the  channel  pair  each  day,  the  numbers  used  are 
shown  in  Table  XII  above.   The  percentage  of  first  run 
messages  and  reruns  that  were  assigned  each  precedence  are 
shown  in  Table  XI  above.   These  percentages  were  the  observed 
aggregate  percentages  for  the  13-16  September  period  being 
simulated. 
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In  the  simulation  the  Broadcast  Control  Station  (BCS) 

handled  all  screen  requests.   The  parameters  used  for  the 

feedback  loop  were: 

i.   Mean  Number  of  Screens  per  Request  =  8 

ii.   Mean  Number  of  Reruns  per  Request  =  1 

iii.   Mean  Screen  Request  Transmission 
Delay  =  120  mins. 

iv.   Mean  Service  Desk  System  Delay  =  75  min. 
To  determine  the  number  of  screens  (messages  to  be  screened) 
that  would  be  handled  each  day  the  regression  equation, 
Y  =  44.36  +  .00065(X),  was  used,  where  X  is  the  number  of 
messages  transmitted  the  previous  day  times  the  number  of 
subscribers.   This  equation  is  taken  from  the  above  described 
regression  analysis  for  FASW  (channel  1)  during  the 
25  August-17  September  period.   Recall  that  during  this 
period  channel  2  was  used  for  rebroadcast  of  channel  1 
traffic.   The  model  determines  the  mean  interorigination 
time  for  screen  requests  by  first  calculating  the  number  of 
screens  and  then  dividing  that  number  by  the  specified  mean 
number  of  screens  per  request.   This  gives  the  number  of 
screen  requests.   By  dividing  the  number  of  screen  requests 
into  1440  minutes  the  result  is  the  daily  mean  interorigina- 
tion time  for  screen  requests. 

To  compare  the  model  results  to  those  actually  observed 
during  the  period  13-16  September  the  model  was  run  in  two 
different  forms.   First  the  model  was  run  with  empirical 
message  length  functions  and  with  empirical  functions  for 
the  service  desk  system  delays  and  the  number  of  reruns  per 
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screen.   Secondly,  the  model  was  run  with  theoretical  func- 
tions for  all  distributions  except  Flash  message  lengths. 
The  Flash  message  length  distribution  was  left  as  an 
empirical  distribution  since  the  size  of  sample  tested  did 
not  allow  us  to  conclude  the  exact  nature  of  its  distribution. 

The  results  of  these  comparison  runs  are  shown  in  Table 
XIII  and  Figure  29  below.   By  comparing  the  average  queue 
lengths  and  the  average  wait  in  queue  we  see  that  the  model 
is  conservative.   We  also  note  that  the  model  with  theoreti- 
cal functions  is  more  conservative  than  the  run  with  empirical 
functions.   In  comparing  the  backlog  curves  in  Figure  29  we 
note  that  the  model  tracks  well.   That  is,  the  model  shows 
peak  backlogs  build  ups  and  declines  at  nearly  the  same  times 
as  were  observed  from  the  actual  data. 

It  should  be  pointed  out  that  we  cannot  say  that  this 
comparison  validates  the  model  since  each  simulation  run  and 
the  actual  observations  for  this  period  are  observations  of 
a  random  process.   Even  if  the  actual  data  could  be  observed 
again  in  the  very  same  circumstances  as  the  period  of  13-16 
September  the  results  would  not  be  the  same,  nor  would 
simulation  model  runs  with  different  seeds  for  the  random 
number  generators  be  the  same. 

F.   MODEL  SIMULATION  RUNS  AND  RESULTS 

1.   Description  of  the  Simulation  Runs 

The  model  with  theoretical  functions  was  used  for 
this  experiment  since  it  was  determined  above  that  it  was 
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more  conservative  than  the  model  with  empirical  functions. 
The  channel  pair  system  was  simulated  under  three  different 
situations,  i)  channel  2  closed  to  first  run  traffic  but 
considered  as  being  used  for  rebroadcast  of  channel  1 
traffic,  ii)  channels  1  and  2  both  transmitting  first  run 
traffic  and  the  rebroadcast  facility  is  not  used,  and  iii) 
channel  2  is  opened  and  closed  during  the  run  by  the  timers 
when  specified  backlogs  were  exceeded.   Additionally  the 
above  situations  were  run  where  reruns  were  allowed  to  go  to 
the  head  of  the  queue  as  opposed  to  going  to  the  back  of 
the  line. 

Eight  simulation  runs  were  made  with  these  config- 
urations.  Four  simulation  runs  were  made  representing  the 
period  of  13-16  September  and  the  other  four  runs  represented 
a  hypothetical  four  days  in  which  the  first  run  traffic  loads 
were  increased  by  200  messages  per  day  over  those  observed 
in  Table  XII  above.   The  runs  made  are  listed  as  follows: 

Run  1.   13-16  SEP.  Channel  2  Closed  for  First  Run 

Run  2.   13-16  SEP.  Channel  2  Open  Whole  Period 

Run  3.   13-16  SEP.  Channel  2  Opened  and  Closed  by  Timers 

Run  4.   Same  as  Run  2  with  Reruns  going  to  the  back  of 

the  line 
Run  5.   High  Loads,  Channel  2  Closed  for  First  Run 
Run  6.   High  Loads,  Channel  2  Open  Whole  Period 
Run  7.   High  Loads,  Channel  2  Opened  and  Closed  by  Timers 
Run  8.   Same  as  Run  7  with  Reruns  going  to  the  Head  of 

the  Line. 

In  runs  1-3,  reruns  were  allowed  to  go  to  the  head 
of  the  queue  and  for  runs  5-7  reruns  were  put  at  the  end  of 
the  queue.   All  of  the  runs  used  the  dirunal  function  shown 
in  Figure  2  8  to  determine  the  hourly  first  run  arrivals.   In 
runs  2-3  and  6-7,  when  channel  2  was  open,  the  number  of 
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screens  were  calculated  by  using  a  different  regression 
equation  and  mean  number  of  screens  per  request  such  that 
the  model  would  simulate  a  higher  number  reruns.   For  runs 
1-4  the  mean  service  desk  system  delay  set  at  75  minutes 
while  for  runs  5-8  it  was  set  at  158  minutes  (see  Table  X) . 

The  allowed  backlog  levels  shown  in  Table  XI  for 
opening  channel  2  in  simulation  runs  4  and  7  were  chosen 
based  on  the  Broadcast  Speed  of  Service  Criteria  shown  in 
Table  I.   By  using  a  familiar  queueing  theory  equation, 
L  =  AW,  where  L  is  the  average  queue  length,  W  is  the  average 
wait  in  queue  and  X  is  message  arrival  rate,  one  can  deter- 
mine the  approximate  backlog  levels  allowed  before  the  Speed 
of  Service  Criteria  are  exceeded  [6,7].   We  note  that  this 
equation  relates  averages  over  a  period  of  time.   Therefore 
the  maximum  allowed  routine  queue  length,  which  is  usually 
the  largest  queue,  was  determined  by  using  this  equation  and 
the  arrival  rate  of  routine  messages  observed  during  the 
13-16  September  period.   The  other  allowed  backlogs  by 
precedence  were  chosen  somewhat  arbitrarily  but  such  that 
the  allowed  backlog  would  not  exceed  on  the  average,  the 
Broadcast  Speed  of  Service  Criteria. 

2 .   The  Results  of  the  Simulation  Runs 

The  results  of  these  eight  simulation  runs  are  shown 
in  Tables  XIV  and  XV  and  Figures  30-32.   In  comparing  the 
average  queue  lengths,  average  wait  in  queue  and  the  maximum 
backlog  observed  plus  the  backlog  curves  the  result  of  the 
simulation  is  obvious.   The  system  configuration  where  both 
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channels  transmit  first  run  traffic  and  the  rebroadcast 
facility  is  not  used  is  superior  to  any  other  configuration. 
Recall  that  the  decision  variable  for  this  problem  is  the 
average  wait  in  queue  of  any  precedence  message  such  that 
this  time  does  not  exceed  the  Broadcast  Speed  of  Service 
Criteria  indicated  in  Table  I.   We  also  note  that  the  primary 
goal  of  the  fleet  broadcast  is  to  make  timely  and  efficient 
delivery  of  traffic  to  all  subscribers.   The  model  shows 
that  there  is  no  question  as  to  the  timeliness  of  traffic 
delivery  when  both  channels  are  running  first  run  traffic. 
So  now  we  must  question  whether  this  configuration  is  also 
the  most  efficient  method.   This  can  be  answered  by  looking 
at  the  feedback  loop  and  the  message  reruns.   Looking  at  the 
results  for  the  high  load  runs  5  and  6  where  channel  2  was 
closed  for  first  run  messages  in  one  run  and  open  for  first 
run  in  the  other  run.   For  these  two  runs  the  mean  time  delay 
in  transmission  of  a  screen  request  was  120  minutes  and  the 
mean  service  desk  system  delay  was  15  8  minutes  thus  on  the 
average  it  took  278  minutes  (4.6  hours)  for  a  screen  request 
to  reach  the  broadcast  position  for  transmission  since  the 
time  it  was  originated.   Now  if  we  chose  the  routine  messages 
for  comparison,  we  note  that  when  channel  2  was  closed  to 
first  run  traffic  the  average  time  it  took  for  a  routine  to 
be  rerun  is  2433  minutes  (40.6  hours)  as  compared  to  326 
minutes  (5.4  hours)  when  channel  2  was  transmitting  first  run 
traffic.   We  also  note  that  15.5%  of  the  total  first  run 
traffic  was  rerun  in  run  6  as  compared  to  only  4.4%  in  run  5. 
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To  fully  explore  this  area  we  further  look  at  the 
comparison  of  the  runs  where  reruns  were  allowed  to  go  to 
the  head  of  the  line.   We  see  that  although  the  reruns  did 
not  have  to  wait  as  long  in  queue  as  they  did  when  they  went 
to  the  back  of  the  line,  we  also  see  that  the  first  run 
traffic  had  a  longer  wait  in  queue  on  the  average  than  when 
reruns  were  put  at  the  back  of  the  line.   This  leads  us  to 
ask  the  question,  "Which  type  of  message  is  most  important?" 
It  would  appear  that  the  first  run  message  is  most  important 
because  it  has  not  been  transmitted  to  any  of  the  subscribers 
copying  the  broadcast  and  hence  it  is  really  undisclosed 
information  while  a  rerun  message  has  most  likely  been 
received  by  the  majority  of  addressees,  if  more  than  one, 
and  hence  it  is  old  information  to  those  subscribers  except 
for  the  subscriber  who  missed  the  message. 
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V.   CONCLUSIONS 

The  conclusions  which  can  be  made  are  that  the  broadcast 
channel  pair  configuration  where  both  channels  transmit  first 
run  traffic  is  the  superior  configuration  for  timely  delivery 
of  traffic.   The  comparison  of  time  delays  for  a  subscriber 
to  receive  a  rerun  show  that  this  configuration  is  not  only 
timely,  it  is  efficient  as  well.   Therefore,  the  answer  to 
the  problem  question  is  that  the  system  should  always  be  run 
with  the  parallel  channels  transmitting  first  run  traffic 
and  that  the  rebroadcast  of  first  run  traffic  on  the  secondary 
channel  should  be  discontinued. 
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APPENDIX  A 
ORIGIN  OF  THE  STUDY 

.  The  problem  on  which  this  study  is  based  was  proposed  as 
a  thesis  topic  to  this  student  in  April  1971  by  Captain  F.  M, 
Snyder,  USN,  then  the  Assistant  Commander  for  Operations  and 
Readiness  to  the  Commander,  Naval  Communications  Command 
(COMNAVCOMM) .   The  study  was  officially  commissioned  as  a 
thesis  study  by  COMNAVCOMM  message  R  021944Z  SEP  71.   A  copy 
of  this  message  is  shown  below  for  reference  and  as  an 
example  of  a  Naval  Communications  teletype  message. 


CZCGAA050 

RTTUZYUW  RHELSAA006  5  2452004-UUU — RUWJAGA 

ZNR  UUUUU 

R  021944Z  SEP  71 

FM  COMNAVCOMM 

TO  RUWJAGA/NAVAL  POSTGRADUATE  SCHOOL 

INFO  RUWNSAA/NAVCOMMSTA  SFRAN 

BT 

UNCLAS    //N00JW// 

FOR  LCDR.  HARLAN  D.  OELMANN 

THESIS  STUDY 

1.  YOUR  THESIS  STUDY  OUTLINE  HAS  BEEN  REVIEWED  AND  ITS 
PURSUIT  IS  STRONGLY  ENCOURAGED.  COMNAVCOMM  REQUESTS 
COPIES  OF  THE  COMPLETED  STUDY  AND  ANY  INTERIM  REPORTS 
THAT  ARE  DEVELOPED. 

2.  DR.  HARDY  (OEG)  WILL  VISIT  NAVCOMMSTA  SANFRAN  DURING 
EARLY  SEPT.,  AND  IT  IS  SUGGESTED  THAT  YOU  MEET  OR  SPEAK 
WITH  HIM  AT  THAT  TIME  IN  ORDER  TO  BECOME  APPRISED  OF 
OEG  EFFORTS  IN  THIS  FIELD. 

3.  SHOULD  COMNAVCOMM  ASSISTANCE  APPEAR  NECESSARY  OR 
DESIRABLE  IN  THE  DEVELOPMENT  OF  THE  STUDY,  REQUEST  YOU 
MAKE  YOUR  NEEDS  KNOWN  VIA  NAVCOMMSTA  SANFRAN. 

BT 
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APPENDIX  B 
PROBLEM  BACKGROUND 

I.   A  BRIEF  EXPLANATION  OF  THE  NAVAL  COMMUNICATIONS  SYSTEM 

A.  ORGANIZATION  OF  THE  NAVAL  COMMUNICATIONS  SYSTEM 

The  Naval  Communications  System  is  composed  of  2  7  Naval 
Communications  Stations  (NAVCOMMSTAs ) ,  4  Naval  Communications 
Units  (NAVCOMMUs)  and  several  Naval  Radio  Stations 
(NAVRADSTAs)  [8].   These  stations  are  located  throughout  the 
world  and  are  organized  into  several  Naval  Communications 
Areas  (NAVCOMMAREAs) .   The  primary  purpose  of  these  stations 
is  to  relay  messages  to  and  from  ships  of  the  fleet  and 
shore  establishments  within  their  area  of  responsibility. 
The  stations  within  a  NAVCOMMAREA  are  organized  under  one 
NAVCOMMSTA  which  is  designated  as  the  Naval  Communications 
Area  Master  Station  (NAVCAMS)  while  the  other  stations  are 
designated  as  Naval  Communications  Area  Local  Stations 
(NAVCALS) .   Although  the  communication  means  which  each  of 
these  stations  use  take  several  forms,  such  as  teletype, 
voice  and  CW,  etc. ,  this  study  is  concerned  only  with  tele- 
type communications  and  the  fleet  broadcast  system.   There- 
fore discussion  will  consider  only  those  areas. 

B.  TELETYPE  COMMUNICATIONS 

It  is  assumed  that  to  most  readers,  teletype  communica- 
tions is  a  familiar  means  of  communication  which  does  not 
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require  detailed  explanation,  except  to  note  that  when 
messages  are  handled  physically  within  a  station  they  are 
handled  in  teletype  tape  or  chad  tape  form.   Hence  in 
further  discussion  when  we  speak  of  messages  being  handled 
in  a  station  it  is  implicitly  meant  that  such  messages  are 
in  tape  form  unless  otherwise  stated. 

1.  World  Wide  Teletype  Communications 

The  Defense  Department  World  Wide  telecommunications 
system  controlled  by  the  Defense  Communication  Agency  pro- 
vides for  long-haul  and  point-to-point  telecommunications 
between  communications  stations  of  the  Armed  Forces  which 
include  the  NAVCOMMSTAs  and  other  stations  within  the 
NAVCOMMAREAs .   This  system  is  known  as  the  AUTODIN  network 
which  is  a  digital  high  speed  message  transmission  system. 

2.  Local  Teletype  Communications 

The  primary  means  of  communications  between  a 
NAVCOMMSTA  and  the  ships  at  sea  is  a  teletype  system  which 
operates  at  the  rate  of  10  0  words-per-minute  (wpm) ,  where 
one  teletype  word  is  equivalent  to  5  teletype  characters. 

There  are  essentially  three  methods  of  teletype 
communication  between  the  ships  at  sea  and  the  shore  based 
NAVCOMMSTA. 

The  first  is  a  ship-shore  telecommunication  link 
which  are  time  shared  circuits  used  by  several  ships  at  sea, 
these  circuits  provide  the  capability  to  send  messages  to  a 
NAVCOMMSTA  for  further  relay  and  delivery. 
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The  second  method  is  also  a  ship-shore  link  which  is 
a  dedicated  circuit,  sometimes  called  a  full  period  termina- 
tion circuit/  between  one  ship  and  the  NAVCOMMSTA.   This 
type  of  circuit  is  normally  used  by  the  larger  ships  such 
as  aircraft  carriers  and  cruisers,  which  normally  have 
staffs  embarked.   The  full  period  termination  circuit  pro- 
vides the  capability  for  these  ships  to  send  messages  to  the 
NAVCOMMSTA  for  further  relay  and  also  to  receive  from  the 
NAVCOMMSTA  messages  which  are  addressed  to  that  particular 
ship  or  embarked  staff. 

The  third  method  is  the  fleet  broadcast  system. 
This  system  is  the  primary  means  of  delivering  traffic  to 
the  fleet  as  a  whole  because  it  provides  for  simultaneous 
delivery  of  messages  to  the  operating  forces  dispersed  over 
a  large  geographical  area.   This  method  also  employs  a 
100  wpm  teletype  system  which  could  be  described  as  a  one 
way  system  because  the  broadcast  system  enables  the  NAVCOMMSTA 
only  to  send  messages  while  the  ships  can  only  receive 
messages . 

The  fleet  broadcast  system  is  composed  of  several 
different  types  of  broadcasts.   These  are  the  fleet  multi- 
channel broadcast,  on  which  this  study  is  based,  a  single 
channel  broadcast  for  ships  that  are  not  equipped  to  copy 
the  multi-channel  broadcast,  the  submarine  broadcast  and  a 
merchant  ship  broadcast.   In  some  NAvCOMMAREAs  these  broad- 
casts are  all  operated  by  one  NAVCOMMSTA  or  Broadcast  Control 
Station  (BCS)  while  in  other  NAVCOMMAREAs  these  separate 
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broadcasts  are  operated  by  two  or  more  NAVCOMMSTAs .   Also 
if  the  capability  of  flexibility  exists,  a  separate  addi- 
tional broadcast  may  be  created  within  a  NAVCOMMAREA  to 
handle  special  traffic  load  situation.   Sometimes  a  NAVCAMS 
may  shift  a  braodcast  to  a  NAVCALS  in  an  emergency  or  to 
maintain  proficiency  among  the  capable  NAVCALS  in  operating 
the  broadcast. 

C.   MESSAGE  PRECEDENCES  AND  SPEED  OF  HANDLING 

Naval  communicaions  messages  are  assigned  one  of  four 
precedences.   Each  message  precedence  has  a  Speed  of  Service 
Criteria  assigned  by  reference  9,  which  is  a  time  standard 
for  a  message  to  take  in  transmission  from  the  origination 
of  the  message  to  the  receipt  of  the  message  by  the  addressee 

These  criteria  are  shown  in  Table  XVI  below.   The  letter 
in  parenthesis  is  the  designator  of  the  precedence. 

TABLE  XVI 
SPEED  OF  SERVICE  CRITERIA 


Precedence 


Average  Speed 
of  Seriice 


Limits 


Flash  (Z) 

Immediate  (0) 
Priority  (P) 
Routine  (R) 


As  fast  as 
humanly  possible 

30  minutes 

3  hours 

6  hours 


Same 

30  minutes-1  hour 

1-6  hours 

3  hours-to  the 
start  of  the  next 
business  day 
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II.   THE  FLEET  MULTI-CHANNEL  BROADCAST  SYSTEM 

There  are  normally  five  fleet  multi-channel  broadcasts 
operated  in  the  Naval  Communications  System.  Each  broadcast 
which  serves  a  particular  ocean  area  is  operated  by  a  desig- 
nated NAVCOMMSTA  or  Broadcast  Control  Station.  These  broad- 
casts have  designations  of  NMUL,  EMUL,  KMUL,  FMUL  and  GMUL , 
where  FMUL  serves  the  Eastern  and  Mid-Pacific  Ocean  area  and 
is  normally  operated  by  NAVCOMMSTA  San  Francisco  [10] . 

A.  CHANNEL  ALIGNMENT  AND  USE 

The  fleet  multi-channel  broadcast  is  a  multiplex  system 
which  has  8  separate  100  wpm  teletype  traffic  channels 
which  are  transmitted  on  a  single  carrier  frequency.   The 
alignment  of  the  traffic  channels  is  such  that  the  overall 
broadcast  traffic  load  to  the  fleet  is  segregated  on  the 
basis  of  ship  type  where  each  subscriber  is  served  by  two 
traffic  channels,  a  primary  channel  and  a  secondary  channel. 
An  example  of  the  overall  broadcast  channel  alignment, 
channel  designations  and  uses  of  each  channel  is  shown  in 
Table  XVII  below. 

Only  channel  pairs  1  and  2,  3  and  4,  and  5  and  7  will  be 
considered  for  this  study. 

B.  THE  RECEPTION  SIDE  OF  THE  BROADCAST 

The  fleet  multi-channel  broadcast  system  for  each  channel 
is  a  guard  system.  This  means  that  each  subscriber  copying  a 
particular  channel  screens  each  sequentially  serial-numbered 
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TABLE  XVII 
FMUL- BROADCAST  CHANNEL  ALIGNMENT 


Channel   Designation 


Use 


1 
2 

3 

4 


FASW 
FSPC 

FALD 

FUSN 

FNSC 

FOP  I 
FHIC 
FMET 


Destroyer  force  primary  traffic  channel 

1  hour  delayed  rebroadcast 'of  channel  1 
or  overload  as  necessary 

Service,  Amphibious  and  Mine  force 
primary  traffic  channel 

1  hour  delayed  rebroadcast  of  channel  3 
or  overload  as  necessary 

Major  warships  (carriers,  cruisers, 
large  destroyers,  including  flagships) 
primary  traffic  channel 

Special  traffic  channel  (Not  of  study 
interest) 

1  hour  delayed  rebroadcast  of  channel  5 
or  overload  as  necessary 

Weather  traffic  channel  (Not  of  study 
interest) 


message  sent  on  that  channel  to  determine  whether  that 
particular  message  is  addressed  to  him  or  not.   If  the  mes- 
sage is  for  the  ship  or  for  the  embarked  staff  the  message  is 
processed.   If  it  is  not  for  that  particular  ship  or  embarked 
staff  the  message  is  ignored  and  filed.   Hence  a  subscriber 
must  perform  two  necessary  actions.  '  He  must  determine  from 
the  message  heading  (which  contains  the  addressees)  whether 
the  ship  or  embarked  staff  is  an  addressee  or  not.   If  the 
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message  is  addressed  to  the  ship  then  it  must  be  copied  with- 
out errors  from  beginning  to  end.   That  is  the  addressee 
must  have  a  clear  enough  copy  that  he  is  able  to  obtain  all 
the  information  transmitted  without  any  doubt  as  to  what 
the  message  says. 

If  a  subscriber  has  missed  part  or  all  of  the  heading 
such  that  he  cannot  determine  whether  the  message  was 
addressed  to  him  or  not,  or  if  he  has  missed  part  or  all  of 
the  text  then  we  will  refer  to  such  a  message  as  a  "missed 
message. " 

A  message  may  be  missed  for  various  reasons,  some  of 
which  are:   environmental  or  atmospheric  conditions,  which 
cause  deterioration  of  the  radio  signal,  equipment  failure 
on  the  ship  or  in  some  cases  error  in  judgment  of  the 
operator  monitoring  the  teletype  printer. 

C.   PROCEDURE  FOR  OBTAINING  MISSED  MESSAGES 

If  a  subscriber  copying  a  particular  channel  has  missed 
one  or  more  messages  he  must  determine  whether  the  message (s) 
was  addressed  to  him  and  if  necessary  he  must  obtain  a 
complete  copy  of  the  message.   He  is  first  encouraged  to 
obtain  a  copy  from  another  ship  copying  the  same  broadcast 
channel,  providing  the  subscriber  has  a  separate  means  of 
communication  with  the  other  ship.   When  the  subscriber 
exhausts  all  local  means  of  obtaining  this  information  his 
next  action  is  to  originate  a  message  which  is  called  a 
"service  or  screen  request  message."   Such  a  message  is 
usually  sent  to  a  NAVCOMMSTA  requesting  them  to  verify  or 
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screen  their  copy  of  the  broadcast  to  determine  whether  the 
message (s)  in  question  is  for  the  subscriber.   The  NAVCOMMSTA 
will  in  turn  reply  by  indicating  which  messages  are  of  no 
concern  to  the  subscriber  and  by  rerunning  on  the  primary 
channel  those  messages  which  were  addressed  to  the  sub- 
scriber.  A  subscriber  which  misses  only  a  portion  of  a 
message  may  service  only  for  that  portion  of  the  message  to 
which  the  NAVCOMMSTA  must  in  turn  originate  a  service 
message  in  reply  giving  the  subscriber  that  portion  of  the 
message  missed. 

It  is  noted  that  in  some  NAVCOMMAREAs ,  the  NAVCOMMSTA 
to  which  a  screen  request  is  sent  may  not  be  the  Broadcast 
Control  Station.   In  these  cases  designated  NAVCOMMSTAs  in 
the  NAVCOMMAREA  also  copy  the  broadcast  and  serve  as 
stations  which  handle  screen  requests  for  a  particular 
channel  of  the  broadcast.   When  these  NAVCOMMSTAs  process  a 
screen  request  their  replys  are  relayed  via  AUTODIN  to  the 
Broadcast  Control  Station  which  in  turn  sends  the  reply  on 
the  appropriate  broadcast  channel.   We  note  that  these  replys 
are  handled  in  the  same  manner  as  first  run  messages. 

D.   HANDLING  OF  SERVICE  MESSAGES  AT  A  NAVCOMMSTA 

Some  NAVCOMMSTAs  will  handle  screen  requests  for  more 
than  one  broadcast  channel  or  in  some  cases  for  more  than 
one  type  of  broadcast  system.   At  each  NAVCOMMSTA  all  such 
broadcast  screen  requests  are  usually  handled  at  one  broad- 
cast service  desk.   The  broadcast  service  desk  handles  these 
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requests  on  the  highest  precedence  first  and  on  a  first  in 
first  out  basis  within  precedences. 

These  service  messages  are  usually  referred  to  as  screen 
requests,  since  the  service  desk  clerk  is  required  to  "screen" 
the  heading  of  the  missed  messages  in  the  printed  copy  form 
of  the  broadcast  channel  log  to  determine  if  the  message  is 
for  the  requesting  subscriber.   Additionally,  to  clarify 
terminology,  each  message  looked  up  is  referred  to  as  a 
"screen."   Thus  a  "screen  request"  may  contain  one  or  more 


"screens . 


After  the  service  desk  clerk  determines  the  result  of  a 
screen  request  he  then  must  make  teletype  tape  copies  of  the 
messages  to  be  rerun,  if  necessary,  which  are  forwarded  to 
the  appropriate  broadcast  channel  position  for  transmission. 
At  some  NAVCOMMSTAs ,  when  the  service  desk  clerk  has  finished 
the  screen  and  before  he  prepares  the  rerun  tapes  he  prepares 
a  short  message  to  the  requester  indicating  which  messages 
are  of  no  concern  to  the  subscriber  and  additionally  states 
that  those  messages  which  are  of  concern  will  be  rerun  later. 
Other  NAVCOMMSTAs  use  the  procedure  where  they  include  this 
information  with  the  first  message  rerun  and  only  after  the 
entire  screen  request  has  been  processed. 

It  is  apparent  that  this  process  is  a  very  time  consuming 
task  and  explains  the  reason  why  subscribers  are  encouraged 
to  seek  other  sources  for  verification  before  asking  the 
NAVCOMMSTA. 
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E.   TRANSMISSION  SIDE  OF  THE  BROADCAST 

Considering  now  the  transmitting  side  of  the  broadcast 
we  look  at  the  NAVCOMMSTAs  internal  handling  of  messages 
which  are  destined  for  transmission  on  one  or  more  of  the 
broadcast  channels.   Messages  arrive  at  the  broadcast  area 
in  teletype  tape  form  and  are,  depending  on  individual 
NAVCOMMSTA  configurations ,  primarily  delivered  to  each  broad- 
cast channel  position  by  two  850  wpm  reperforators  called 
"burpees."   The  traffic  delivered  by  the  burpees  is  traffic 
which  has  been  received  by  the  NAVCOMMSTA  from  the  AUTODIN 
network  or  from  100  wpm  teletype  circuits  which  terminate  at 
the  NAVCOMMSTA  from  any  number  of  sources.   The  burpees  are 
fed  from  these  various  sources  through  a  small  computer 
system  called  the  Multiple  Address  Processing  Unit  (MAPU) . 
This  system  reads  the  header  of  a  message  which  contains 
address  codes  as  to  where  a  message  is  destined  within  the 
NAVCOMMSTA  and  routes  that  message  to  the  appropriate  relay 
positions  by  delivering  a  teletype  tape  through  the  burpees. 
In  addition  to  the  burpee  deliveries  some  NAVCOMMSTAs  hand 
deliver  screen  request  answers,  reruns  and  other  messages  to 
the  broadcast  position.   The  screen  request  answers  and 
reruns  are  then  usually  put  at  the  head  of  the  line  or  queue 
for  transmission. 

Messages  delivered  to  the  broadcast  position  are  put  in 
a  grid  (queue)  by  precedence  and  in  order  of  arrival  within 
each  precedence.  The  messages  are  then  taken  from  the  grid 
and  put  in  the  channel  transmitter  distributor  (TD) ,  a 
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device  which  reads  the  teletype  tape  for  transmission.   The 
order  in  which  the  messages  are  taken  from  the  grid  for 
transmission  is  on  the  highest  precedence  waiting  and  first 
in  first  out  (FIFO)  within  a  precedence  basis.   This  process 
at  some  NAVCOMMSTAs  is  done  by  one  man  physically  performing 
the  queueing  process  while  another  man  pulls  the  tapes  from 
the  grid  and  places  them  into  the  TDs  for  transmission.   At 
other  NAVCOMMSTAs,  properly  equipped,  this  entire  process  is 
handled  by  a  small  computer  system,  known  as  the  "tight 
tape"  system.   In  either  case  the  queueing  and  transmission 
flow  is  the  same  for  the  purpose  of  this  study. 
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APPENDIX  C 
GOODNESS  OF  FIT  DATA  ANALYSIS  COMPUTER  PROGRAM 

I.   DESCRIPTION  OF  CAPABILITY 

The  Goodness  of  Fit  Data  Analysis  Program  accepts  any 
number  of  data  samples  of  varying  sizes  and  performs 
prescribed  calculations  and  produces  various  graphs  and 
tables  which  are  used  to  analyze  the  fitting  of  an  empirical 
distribution  with  a  selected  theoretical  distribution.   The 
program  is  written  to  fit  the  empirical  data  to  two  different 
distributions.   These  two  distributions  are  the  Exponential 
and  the  Hyperexponential  probability  distributions  in  either 
their  regular  form  or  the  right  shifted  form. 

The  main  program  uses  five  different  subroutines  which 
are  named  SORT,  STAT,  HISTO,  PLOTP  and  UTPLOT,  the  last  two 
subroutines  being  used  in  conjunction  with  each  other.   The 
main  program  with  the  subroutines  first  sorts  the  data  in 
ascending  order  and  prints  out  the  data  in  this  form.   Then 
the  program  calculates  statistical  parameters  and  prints  out 
the  parameter  values  which  are  later  used  in  subsequent 
calculations.   The  next  feature  provides  for  sorting  the 
data  into  frequency  classes  and  prints  a  histogram  of  the 
data  according  to  the  number  and  widths  of  intervals  as 
prescribed  by  the  programmer.   The  program  then  calculates 
and  plots  the  log  of  the  tail  distribution  for  both  the 
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empirical  and  theoretical  distributions  being  tested.   A 
description  of  these  calculations  along  with  other  mathe- 
matical descriptions  are  more  fully  explained  below.   The 
respective  Cumulative  Distribution  Function  (CDF)  values  are 
calculated  and  the  absolute  difference  of  the  distribution 
points  is  determined  and  printed.   The  absolute  difference 
provides  the  user  this  information  for  use  in  performing 
the  Kolmogorov-Smirnov  (K-S)  Goodness  of  Fit  Test  on  the 
data.   The  cumulative  distribution  functions  are  plotted  to 
facilitate  a  graphical  interpretation  of  the  fit  of  the 
distributions.   Finally  the  data  is  analyzed  by  performing 
the  calculation  of  the  Chi-Square  Statistic  and  printing  out 
a  table  of  results  for  use  in  performing  the  Chi-Square 
Goodness  of  Fit  Test. 

Features  of  the  program  provide  a  flexibility  of 
selecting  only  certain  calculations  to  be  performed  and/or 
certain  outputs  to  be  printed.   These  features  of  selection 
are  fully  documented  in  the  computer  program  shown  in  this 
Appendix. 

This  Appendix  is  not  meant  to  be  a  user's  guide  for  the 
program  but  rather  to  explain  and  reference  the  mathematical 
calculations  used  in  the  analysis  program. 

II.   THE  THEORETICAL  DISTRIBUTIONS 

A.   THE  EXPONENTIAL  DISTRIBUTION 

The  exponential  density  function/  cumulative  distribution 
function  and  the  associated  parameters  are  given  below. 
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f  (x)  = 


Xe 
0 


-X  (x-a) 


=  0 


,  x>_a,  where  a  is  the  right 
shift  parameter. 

,x<a, 


F(x)  =  i-e~X(x"a)    ,x>a, 


,x<a. 


E(X)  =  ±  +  a, 


y-a' 


=  VAR(X)  = 


2  " 


B.   THE  HYPEREXPONENTIAL  DISTRIBUTION 

The  hyperexponential  distribution  shown  below  is  a 
mixture  of  two  exponential  distributions.   The  observations 
originate  in  two  populations  only,  with  probabilities  p  and 
1-p  and  parameters  2Xp  and  2X (1-p) ,  respectively  [11]. 

The  hyperexponential  density  function,  cumulative  dis- 
tribution function  and  the  associated  parameters  are  given 
below. 

£  ,    .         -.2,  -2pX  (x-a)  ,  ~/n   .,  -2  (1-p)  X  (x-a) 

f(x)  =  2p  Xe   ^      +  2(l-p)Xe     trt  ,      x>a, 

=  0  ,   x<a, 


F(x)  =  p 
=  0 


'1_e-2pX(x-a)l+  (1_p) 


'] 


1_e-2(l-p)X(x-a) 


,    x>a, 
,  x<a, 


where  a  is  the  right  shift  parameter. 


y   =  E(X)  =  ±  +  a, 


X   = 


V-a' 
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a  = 


P  = 


■\?n d  ^y^     — 

] 

l 

1              1 

X2 

_2p(l 
2 

-P) 

h 

2         2 

1       a2 

y 

+    1- 

III.   STATISTICAL  AND  GRAPH  CALCULATIONS 

A.   SAMPLE  STATISTICAL  PARAMETERS 

The  analysis  program  has  a  feature  in  which  preselected 
values  for  the  mean,  the  right  shift  parameter  a  and  the 
hyperexponential  parameter  p  can  be  read  in  instead  of  having 
the  program  calculate  these  estimates  from  the  data.   When 
this  feature  is  not  used  some  of  these  parameters  are  cal- 
culated by  the  subroutine  STAT  using  the  familiar  formulas 
shown  below. 


N 


N 


XBAR  =  X 


=  Iv 


VAREST 


i=l 

N 


=  N^T  1      (Xi"  *>  ' 
i=l 


We  note  that 


XBAR  =  V    =  i  +  a  =  ^^   +  A, 


hence , 


LAMEST  = 


(XBAR  -  A) 


The  capitalized  names  are  those  used  in  the  program  for 
these  various  parameters. 
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The  estimate  of  the  hyperexponential  parameter  p  is 
calculated  using  the  same  formula  shown  above  except  that 
XBAR  =  y  and  VAREST  =  a  . 

The  right  shift  parameter  a  can  also  be  specified  to  be 
the  minimum  order  statistic  of  the  data  sample. 


B.   A  GRAPHICAL  PROCEDURE  OF  TESTING  FOR  EXPONENTIALITY 

When  comparing  empirical  data  to  the  exponential  distri- 
bution, for  goodness  of  fit,  this  graphical  technique 
provides  an  appealing  method  of  either  supporting  or  denying 
the  more  purely  statistical  comparisons  to  be  described 
below. 

If  the  empirical  data  is  really  distributed  exponentially 
with  a  cumulative  distribution  function  as  shown  above  then 
the  following  graphical  comparison  can  be  made  [12],   For 
the  exponential  distribution  the  log  of  the  reciprocal  of 
the  tail  distribution  (called  the  log  of  the  tail  for 
brevity)  is: 

,       1      (x-a)         .      ^  „ 
y  =  log  j_F(x)  =  — jj-S   assuming  y  >  0. 

Thus  if  we  plot  y  against  x,  we  get  a  straight  line  with  a 
slope  of  —  and  the  line  will  cut  the  x-axis  at  the  point  x=a. 

If  we  have  n  data  points  in  the  empirical  data,  we  can 
sort  the  n  points  in  ascending  order  forming  the  order 

statistics,  i.e.,  x,  <_   x?  <_ £  x  .   The  empirical 

cumulative  distribution  function  is  defined  as 

F(x.)  =  -  . 
l    n 
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If  we  let  y  =  log  i_F/x  \  •   Then  by  plotting  y.  against  x. 
the  plotted  points  can  be  compared  to  the  fitted  straight 
line  as  a  test  of  the  exponential  assumption. 

It  is  noted  that  the  program  has  the  option  of  using 
a  preselected  mean  or  using  the  sample  mean  estimate  in 
calculating  the  straight  line.   Thus,  this  graphical  com- 
parison can  be  used  to  make  an  independent  estimate  of  the 
mean  by  determining  the  slope  of  the  empirical  data.   Then 
the  data  can  be  rerun  with  a  preselected  mean  to  attempt 
a  better  fit. 

When  the  log  of  the  tail  distribution  is  used  with  the 
hyperexponential  comparison  we  get  a  decreasing  failure  rate 
function  for  the  hyperexponential  distribution.   That  is  we 
get  a  concave  curve.   In  this  case  the  graphical  comparison 
can  also  be  used  to  determine  how  close  the  two  compare. 

This  method  of  distribution  fitting  can  be  used  when 
empirical  data  is  analyzed  for  the  first  time.   Through  the 
use  of  the  log  of  the  tail  graph  we  can  determine  how  the 
distribution  behaves  and  find  out  whether  it  has  a  decreas- 
ing or  increasing  failure  rate  and  thus  plan  further  analysis 
on  the  data. 

C.   THE  KOLMOGOROV-SMIRNOV  GOODNESS  OF  FIT  TEST 

The  Kolmogorov-Smirnov  Goodness  of  Fit  Test  is  a  familiar 

test  and  therefore  will  not  be  described  here  [13] . 

It  is  noted  that  the  program  only  calculates  the  absolute 

difference  between  the  F(x)  values  of  the  two  distributions 

and  prints  the  results.   Therefore  the  user  must  scan  the 
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points  to  find  the  maximum  difference  and  then  enter  the 
appropriate  Kolmogorov-Smirnov  table  to  determine  the 
results  of  this  test. 

The  Kolmogorov-Smirnov  test  is  regarded  by  most 
references  as  a  more  powerful  test  than  the  Chi-Square 
Goodness  of  Fit  Test  but  it  has  one  drawback,  that  being 
that  the  Kolmogorov-Smirnov  statistic  becomes  inflated  when 
the  data  has  large  groups  of  points  which  are  of  equal  value. 

D.   THE  CHI-SQUARE  GOODNESS  OF  FIT  TEST 

The  Chi-Square  Goodness  of  Fit  Test  is  also  a  familiar 
test  and  will  be  described  only  briefly  [12,14,15]. 

It  is  well  known  that  the  chi-square  goodness  of  fit 
test  has  several  drawbacks.   Among  them  are  its  large  sample 
character  and  the  dependence  upon  the  careful  choice  of  the 
number  and  the  widths  of  the  intervals.   Therefore  this  test 
should  be  used  with  caution.   For  further  discussion  of 
these  characteristics  see  reference  14. 

It  is  also  noted  that  the  program  has  the  capability  of 
either  reading  in  a  specified  mean  or  estimating  it.   This 
will  affect  the  number  of  degrees  of  freedom  used  when 
testing  the  chi-square  statistic  against  the  appropriate 
table  values.   The  chi-square  goodness  of  fit  test  with  an 
estimated  mean  uses  (k-2)  degrees  of  freedom  or  (k-1)  degrees 
of  freedom  when  the  mean  is  preselected,  where  k  is  the 
number  of  intervals  (or  classes)  used  in  the  test,  and  where 
the  x-axis  is  divided  into  the  intervals  x,  <_   x?  < — xk-l 
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The  chi-square  statistic  was  calculated  in  the  following 
manner.   Let  o.  be  the  observed  number  of  data  points  in  the 
i    interval  and  let  e.  be  the  expected  number  of  data  points 
in  the  i    interval,  where  e.  =  np .  and  where  n  is  the 
number  of  data  points  in  the  sample  being  tested  and  p.  is 
calculated  as  indicated  below  for  the  appropriate  distribu- 
tion. 

For  the  exponential  distribution: 
x.  +  l 


i=/\e-<t 


P^       = 


-X(x.-a)      -X(x.J^1-a) 

-a)  , .  1  l+l 

dt  =   e  -e 


1^  X    /    •     •     •     f  jC~"  rib    f 


and 


.-/ 


,    -X(t-a),,.  -x^xk_i-a)         . 

Xe  dt  =   e  ,      l   =   k-1. 


^k-1 


For  the  hyperexponential  distribution: 

xi+l 

"V,    2,    -2pX(t-aK     0/1       x,    -2  (1-p)  X  (t-a) 
2p    Ae      ^  +    2(l-p)Xe  c 


P,-    = 


X. 

1 


dt 


-.[. 


-2pX(xi-a)       -2pX(xi+1-a) 


-e 


-2  (1-p) X (Xi-a)      -2 (1-p) X (xi+1-a) 


+    (1-p)    e 


and 


i 


=   pe 


2p2Xe"2PX  (t"a)  +  2  (1-p)  Xe"2  (1-P»  X  {t'^ 

_ 

k-1 
-2pX (xk_x-a) 


dt 


-2(l-p)X(x        -a) 
+    (1-p)  K"X         ,       i=k-l. 
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The 


value  e.  was  restricted  to  the  condition  of  having 


e.  >  5  for  each  interval  [14].   If  this  condition  is  not 

1  — 

satisfied  then  the  observations  of  that  interval  are  added 
to  the  previous  interval  and  the  appropriate  integral  is 
recalculated.   This  is  repeated  until  this  condition  is 
satisfied. 

Finally  the  chi-square  statistic  is  determined  by  this 
formula 

2       £   (o,-e,)2 


stat 

1=1 

This  statistic  is  then  compared  to  the  appropriate 
chi-square  table  value  to  determine  the  results  of  the  test 
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IV.   DATA  ANALYSIS  COMPUTER  PROGRAM 


C 

C  GOODNESS  OF  FIT 

C  DATA  ANALYSIS  PROGRAM 

C 

c 

C   PROGRAM  DISCRIPTION: 

C 

C   SORTS  DATA  IN  ASCENDING  ORDER  AND  PRINTS  DATA 

C   CALCULATES  PARAMETERS:  MEAN,  VARt  STD.DEV.,  LAMEST 

C   AND  HYPEREXPONENTIAL  PROPORTION  PARAMETERS  PI  AND  P2 

C 

C   HISTOGRAMS  SORTED  DATA 

C 

C   CALCULATES  AND  PLOTS  THE  NATURAL  LOG  OF  THE  RECIPROCAL 

C   OF  THE  TAIL  DISTRIBUTION  FOR  BOTH  THE  EMPIRICAL  AND 

C   THE  THEORETICAL  DISTRIBUTIONS 

C 

C       CALCULATES    THE    C.D.F.'S    FOR    BOTH    THE    EMPIRICAL    AND 

C       THE    THEORETICAL    DISTRIBUTIONS    AND    DETERMINES    THE 

C       ABSOLUTE    DIFFERENCE    BETWEEN    THE    TWO    C.D.F.'S    FOR    USE 

C       IN    THE    KOLMOGOROV-SMIRNOV       GOOONES    OF    FIT    TEST 

C 

C   PLOTS  THE  C.D.F.'S  OF  THE  EMPIRICAL  AND  THEORETICAL 

C   DISTRIBUTIONS  ON  THE  SAME  GRAPH 

C 

C   CALCULATES  CHI-SQUARE  STATISTIC 

C   PRINTS  A  TABLE  OF  THE  CHI-SQUARE  CALCULATION  RESULTS 

C 

REAL  X(450),Y(450,4) ,Z(450) 

REAL  YY(450) ,XR(450) 

REAL  FREQ(IOO) ,RANGE( 100) 

REAL  FREQ1 (20) , R ANGE1 ( 20 ) , T ITL E ( 20 ) 

REAL  P(100),E(100) ,XX( 100,3) 

REAL  INC, LAMEST 

INTEGER  TYDIST,RS 

INTEGER  DPS,PPS,HPS,LPS,KPS,CPS,CSPS 
C 

C   DEFINE:   NUMTIM — NR  OF  SETS  OF  DATA  TO  BE  ANALYZED 
C 
C 

C   DEFINE:   TYPE  OF  THEORETICAL  DISTRIBUTION  TO  BE  USED 
C  TO  FIT  DATA 

C 

C  TYDIST— FOR  EXPONENTIAL  DIST  SET  TYDIST=0 

C  FOR  HYPEREXPONENTIAL  DIST  SET  TYDIST=1 

C 

C   DEFINE:   RIGHT  SHIFT  OPTION 
C 

C  RS — REGULAR  OR  RIGHT  SHIFTED  DIST  WITH  DEFINED 

C  VALUE  FOR  AA,  I.E.,  REGULAR  DIST  AA=0.0 

C  SET  RS=0  OR  READ  SPECIFIED  VALUE  FOR  AA  BELOW 

C  RIGHT  SHIFT  SET  RS=1  THEN  AA=MIN  ORDER  STAT 

C 

READ(5,1001)NUMTIM,TYDIST,RS 
1001  F0RMAT(3I5) 
C 

C   DEFINE:   INI NR  OF  INTERVALS  USED  IN  (FIRST)  HISTOGRAM 

C  AND  FOR  INITIAL  CHI-SQUARE  CALCULATION 

C  (MAX=20) 

C  IN IN  IS  THE  TOTAL  NR  OF  INTERVALS  (MAX=40) 

C  TO  BE  USED 

C  IF  IN  GT  20  THIS  PROVIDES  TWO  HISTOGRAMS 

C  IF  ONLY  CNE  HISTOGRAM  DESIRED  SET  IN=IN1 
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C  INCD — INCREMENT    DETERM INATOR .       USED    TO    CONTROL 

C  THE    WIDTH    OE    EACH    HISTOGRAM    AND    CHI-SQUARE 

C  INTERVAL,     I.E.,     I  NC=  (B-A)  /  I  NICD 

C  A,B LOW    AND    HIGH    RANGES    OF    HISTOGRAM(S) 

C  INCS — IF    GT    40     INTERVALS    ARE    DESIRED    FOR 

C  CHI-SQUAPE    CALCULATION     INDICATE    NUMBER 

C  DESIRED    (MAX=100). 

C  IF  INCS  NOT  DEFINED  SET  INCS=0 

C  INCSD-INCREMENT  DETERMINATOR  FOR  CHI-SQUARE 

C  INTERVAL  WIDTH.  NOTE:  MUST  BE  DEFINED  WHEN 

C  INCS, C, AND  D  ARE  DEFINED. 

C  CD LOW  AND  hIGH  RANGES  FOR  CHI-SQUARE 

C  INTERVALS  GT  40 

C  IF  INCS  NE  0  THEN  C  AND  D  MUST  BE  DEFINED 

C 

C  READ  STATEMENTS  1101  AND  1201  CAN  BE  MOVED  INSIDE  OF 

C  9999  DO  LOOP  TO  ALLOW  CONTROL  CHANGES  ON  INDIVIDUAL  DATA 

C  DECKS 

C 

READ(5,1101)IN1,IN,INCD,A,B,INCS,INCSD,C,D 
1101  FORM AT (3  1 5, 2F 1 0.2, 21 5, 2F 10. 2) 
C 

C  DEFINE:  PRINT  SUPPRESSION  OPTION.   IF  PRINTOUT  OF 

C  INDICATED  INFORMATION  IS  DESIRED  SET=0,  IF  NOT 

C  DESIRED  SET=1. 
C 

C  DPS — PRINT  OF  SORTED  DATA 

C  PPS — PRINT  OF  DATA  PARAMETERS 

C  HPS — PRINT  OF  HISTOGRAM(S) 

C  LSP — PLOT  OF  LN  OF  TAIL  DIST. 

C  KPS — PRINT  OF  CDF  VALUES  AND  KLOMCGOROV-SMI RNOV 

C  STATISTIC 

C  CPS — PLOT  OF  C.D.F.'S 

C  CSPS-PRINT  OF  CHI-SQUARE  STATISTIC 
C 

READ(5,1201)DPS,PPS,HPS,LPS,KPS,CPS,CSPS 
1201  F0RMAT(7I5) 
C 

C  START  DO  LOOP  WHICH  PROCESSES  EACH  SET  OF  DATA 
C 

DO  9999  L=1,NUMTIM 
C 

C  READ  IN  TITLE  OF  DATA  SET 
C 

READ(5,2001)TITLE 
2001  FORMAT (20A4) 
C 

C  DEFINE:  SAMPLE  SIZE  N 
C 

C  DEFINE:  OPTION  ALLOWS  READING  IN  MEAN  VALUE  TO  BE  USED  IN 

C  CDF,  LN  TAIL  DIST  AND  CHI-SQUARE  CALCULATIONS 

C  OPTION  ALSO  ALLOWS  SPECIFYING  VALUES  FOR  AA 

C  OR  FOR  SPECIFYING  VALUES  FOR  PI  AND  P2  IN  HYPER- 

C  EXPONENETIAL  DISTRIBUTION 
C 

C  RXBAR — READ  IN  MEAN  VALUE 

C  AA — OFFSET  SPECIFIED  FOR  EITHER  REGULAR  (AA=0.0) 

C  DIST  OR  RIGHT  SHIFTED  DIST  AA=READ  IN  VALUE 

C  PI, P2— SPECIFY  HYPEREXP  PARAMETERS   NOTE:  P2=1-P1 

C  IF  NOT  USED  SET  =0 
C 

READ(5,3001)N,RXBAR,AA,P1,P2 
3001  F0RMAT(I5,F10.2,F5.2,2F5.4) 

C   READ  IN  DATA  SET 
C 

READ(5,4001) (X(I  ),  1=1, N) 
4001  FORMAT(16F5.0) 
C 
C   LOAD  RANGE  VECTOR 


C 


INC=(B-A)/INCD 


124 


RANGE(1)=A+INC 
DO  1000  1=2, IN 
RANGE  CI )=RANGE(I-1)+INC 
1000  CONTINUE 
C 

C   LOAD  A  LARGE  NUMBER  IN  RANGE(IN)  SO  WHEN  PRINTED  OUT  IN 
C   F5.0  FORMAT  OR  SIMILAR  FORMAT  IT  WILL  PRINT  AS  *****  TO 
C   INDICATE  INFINITY 
C 

RANGE( IN)=10000.0 
C 

C   SORT  DATA  IN  ASCENDING  ORDER 
C 

CALL  SORT(X,N) 
C 

C   CALCULATE  PARAMETERS 
C 

CALL  STATU, N,XBAR,VAR,S) 
C 

C   AA  DEFINES  THE  OFFSET  IS  A  RIGHT  SHIFTED  DIST  IS  USED 
C 

IF(RS.EQ.O)GO  TO  1310 
AA=X(1 ) 
1310  LAMEST=1.0/(XBAR-AA) 
RLAM=0.0 

IF(RXBAR.EQ.0.0)G0  TO  1610 
RLAM=1.0/(RXBAR-AA) 
1610  IF(TYDIST.EQ.O)GO  TO  1710 
IF(P1.NE.0)G0  TO  1710 
C 

C   CALCULATE  PI  AND  P2  PARAMETERS  FOR  HYPEREXPONENTI AL  DIST 
C   IF  REQUIRED 
C 

P1=0.5-0.5*SQRT( 1.0-2.0/(VAR/XBAR**2+1.0) ) 
P2=1.0-P1 
1710  IF(DPS.EQ.1)G0  TO  1810 
C 

C   WRITE  SORTED  DATA  AND  PARAMETERS 
C 

WRITE(6,1002)TITLE,i  X( I) ,I=1,N) 
1002  FORMATC  1>  ,T10, 'DATA  POINTS  AND  ESTIMATED  PARAMETERS' 

X*  FOR: ■ ,20A4//(10F11.2) } 
1810  IFCPPS.EQ.  DGO  TO  1910 

IF(DPS.EQ.1)WRITE(6,2002)TITLE 
2002  FORMATt ■ 1'  ,T10, ■ ESTIMATED  PARAMETERS  F0R:',20A4) 

WRITE(6,2102)XBAR, VAR, S, LAMEST ,N , RXBAR ,RLAM, PI , P2 ,AA 
2102  F0RMAT(//T7,,XBAR=',F11.2, '    VAREST= • , F 15 . 2, 
X'    STD.DEV.=' ,F11.2, 

X'    LAMEST=«,F11.7, •    SAMPLE  SI  ZE  =  '  ,  I  4/T7,  •  RXBAR=  •  , 
XF11.2,'    RLAM=f ,F11.7, •    P1=',F6.4,'    P2=',F6.4, 
X«    AA=',F6.2) 
1910  IF(RXBAR.EQ.0.0)G0  TO  2010 
XBAR=RXBAR 
LAMEST=RLAM 
C 

C   ZERO  OUT  FREQ  VECTOR 
C 
2010  DO  2000  1=1, IN 

FREQU  )  =  0.0 
2000  CONTINUE 
C 

C   DETERMINE  FREQ  OF  DATA  IN  INTERVALS  FOR  HISTOGRAM  AND 
C   CHI-SQUARE  CALCULATION 
C 

Nl  =  IN-2' 
DO  3000  J=1,N 

IF(X(J).LE.RANGE(1))FREQ(1)=FREQ(1)+1 
IF(X(J).GT.RANGE(IN-1))FREQ(IN)=FREQ(IN)+1 
DO  3100  1=1, Nl 

IF(X(  J)  .GT.RANGEd  ).AND.X(  J)  .LE.RANGE(  1  +  1)  ) 
XFREQ(I+1)  =  FREQ(  I+D+l 
3100  CONTINUE 
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3000    CONTINUE 

IF(HPS.EQ.1)G0    TO    3900 
C 

C       PRINT    OUT    FIRST    HISTOGRAM 
C 

CALL    HISTOdNl,FREQ, RANGE, TITLE) 
C 

C       DETERMINE     IF    MORE    THAN    ONE    HISTOGRAM    IS    DESIRED 
C 

IF(  IN.EQ.INDGO    TO    3900 

IN2=IN-IN1 

NSTART=IN1+1 

DO  3200  J=NSTART,IN 

RANG El (J-IN1)=RANGE(J) 

FREQH  J-IN1)=FREQ(  J) 
3200  CONTINUE 

CALL  H ISTO U N2, F RE QltRANGEl, TITLE) 
C 

C   CALCULATE  VALUES  FOR  PLOT  OF  NATURAL  LOG  OF  THE 
C   RECIPROCAL  OF  THE  TAIL  DISTRIBUTION 
C 

C   CALCULATE  LN  OF  TAIL  DIST  FOR  EMP  DIST   PUT  IN  Y(I,1) 
C   CALCULATE  EXPONENTIAL  F(X)  AND  PUT  IN  Yd, 2) 
C   CALCULATE  STRAIGHT  LINE  BENCHMARK  FOR  GRAPH  PUT  IN  Y(I,3) 
C   NOTE:  CALCULATION  OF  THE  STRAIGHT  LINE  IS  NECESSARY  ONLY 
C   WHEN  TESTING  A  THEOR.  DISTRIBUTION  OTHER  THAN  THE   FXP 
C   CALCULATE  LN  OF  TAIL  DIST  FOR  THEOR  DIST  PUT  IN  Yd, 4) 
C   CALCULATION  MADE  FOR  N-l  DATA  PTS  TO  AVOID  TAKING  LNIO) 
C 
3900  NS=N-1 

ARG1=N 

IF(TYDIST.EQ.1)G0  TO  4010 

DO  4000  1=1, NS 

ARG2=N-I 

ARG=ARG1/ARG2 

Y( I, 1 J=ALOG(ARG) 

Y(  I,2)=1.0-EXP(-(LAMEST*(X(I )-AA))  ) 

Y( I,3)=(X(I)-AA)/XBAR 

Y( I ,4) =AL0G(1.0/(1.0-Y(I,2)  )  ) 
4000  CONTINUE 
C 

C   INCLUDE  LAST  DATA  POINT  IN  THEOR  CDF  FOR  CDF  PLOT  BELOW 
C 

Y( N, 2 )  =  1.0-EXP(-( LAMEST* (X(N)-AA)  J  ) 

GO  TO  4110 
4010  DO  4050  1=1, NS 

ARG2=N-I 

ARG=ARG1/ARG2 

Yd,  1)=AL0G(ARG) 

Y( I,2)=P1*(1.0-EXP(-2.0*P1*LAMEST*(X(I  )-AA)  )  )  + 
XP2*( 1.0-EXP(-2.0*P2*LAMEST*(X(I )-AA) )) 

Y( I,3)=(X( I )-AA)/XBAR 

Y(I,4)=AL0G(1.0/(1.0-Y(I,2)  )) 
4050  CONTINUE 
C 

C   INCLUDE  LAST  DATA  POINT  IN  THEOR  CDF  FOR  CDF  PLOT  BELOW 
C 

Y(N,2)=Pl*d.0-EXP  (-2.0*P1*LAMEST*(X(I  )-AA)  )  )  + 
XP2*(1.0-EXP(-2.0*P2*LAMEST*(X( I  )-AA)  )) 
C 

C   PLOT  LN  OF  TAIL  DIST 
C 
4110  IF(LPS.EQ.1)G0  TO  5110 

WRITE(6,1102)TITLE 
1102  FORMAT {• i« ,T10, ■ PLOT  OF  LN ( 1 .0/ 1 .0-F ( X ) )  FOR:   ', 
X20A4//) 

DO  4100  1=1, NS 

YY( I)=Y(I,3) 
4100  CONTINUE 

CALL  PL0TP(X,YY,NS,1) 

DO  4200  1=1, NS 

YY( I )=Y(I ,1) 
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4200  CONTINUE 

CALL  PL0TP(X,YY,NS,2) 
DO  4300  1=1, NS 
YY( I )=Y( 1,4) 
4300  CONTINUE 

CALL  PL0TP(X,YY,NS,3) 
C 

C   CALCULATION  OF  THE  EMPIRICAL  C.D.F.  AND  THE  STATISTICS 
C   FOR  THE  KOLMOGOROV-SMIRNOV  TEST 
C   WRITE  C.D.F.  VALUES  AND  THE  K-S  STATISTICS 
C 

DEM  =  N 

DO  5000  1=1, N 
Y( I,1)=I/DEM 
Z( I)=A3S(  YU ,1)-Y( 1,2) ) 
5000  CONTINUE 
5110  IF(KPS.EQ.1)G0  TO  6110 

WRITE(6, 3002)TITLE 
3002  FORMAT (• 1» ,T10, ■ C.D.F.  VALUES  AND  KOLMOGOROV-SMIRNOV  ' 
X'STATISTIC  FOR: ' //T10 , 20A4//T13 , «X« ,T20, «EMP  DIST   « 
X'EST  DIST    K-S  STAT'/) 
DO  6000  1=1, N 

WRITEC 6,4002 )X( I  ) , (Y( I  ,J)  ,J  =  1,2)  ,Z(I ) 
4002  F0RMAT(T9,F7.2,T20,F8.4,2F11.4) 
6000  CONTINUE 

6110  IFtCPS.EQ.UGO  TO  9110 
C 

C   GRAPH  THE  EMP  AND  THEORETICAL  C.D.F.S 
C 

WRITE(6,4202)TITLE 
4202  FORMAT!' 1'  ,T10,  '  EMPIRICALU  )  AND  THEORET  IC  AL  (  +  )  CDF  • 
X'FOR: « ,20A4//) 
C   ADD  0.0  IN  X  AMD  YY  VECTORS  SUCH  THAT  THE  CDF  PLOT  WILL 
C   AUTOSCALE  FROM  0.0 
N2=N+1 
XR(1 )=0.0 
YY(1)=0.0 
DO  7000  1=2, N2 
XR( I )=X( 1-1) 
YY(I )=Y( 1-1,1) 
7000  CONTINUE 

CALL  PL0TP(XR,YY,N2, 1) 
DO  8000  1=2, N2 
YY(I )=Y( 1-1,2) 
8000  CONTINUE 

CALL  PL0TP(XR,YY,N2,3) 
C 

9110  IF(CSPS.EQ.1)G0  TO  9999 
C 

C   REDEFINE  AND  RELOAD  RANGE  AND  FREQ  VECTORS  FOR  CHI-SQUARE 
C   CALCULATION.   IF  INCS=0  THEN  BYPASS  ROUTINE 
C 

IF( INCS.EQ.O)GO  TO  0099 
C 

C   STORE  VALUE  OF  IN  FOR  NEXT  DECK  OF  DATA  TO  BE  PROCESSED 
C 

TEMPIN=IN 
IN=INCS 
C 

C   LOAD  RANGE  VECTOR 
C 

INC=(D-C)/INCSD 
RANGE( 1)=C+INC 
DO  9000  1=2, IN 
RANGE(  I  )=RANGE(  I-D  +  INC 
9000  CONTINUE 
C 

C   LOAD  A  LARGE  NUMBER  IN  RANGE(IN)  SO  WHEN  PRINTED  OUT  IN 
C   F5.0  FORMAT  OR  SIMILAR  FORMAT  IT  WILL  PRINT  AS  *****    TO 
C   INDICATE  INFINITY 
C 

RANGEt  IN)  =  10000.0 
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c 

C   ZERO  OUT  FREQ  VECTOR 
C 

DO  9100  1=1, IN 

FREQU  )=0.0 
9100  CONTINUE 
C 

C   DETERMINE  FREQ  OF  DATA  IN  INTERVALS  FOR  CHI-SQUARE 
C   CALCULATION 
C 

Nl-IN-2 

DO  9200  J=1,N 

IF(X(J).LE.RANGE(1 )  )  FREQ( 1 ) =FREQ ( 1 )  +  l 

IF(X(J) .GT.RANGE( IN-1) )FREQ{ IN )=FREQ ( I N ) +1 

DO  9300  1=1, Nl 

IF(X( Ji.GT.RANGE (I ).AND.X( J) .LE. RANGE ( 1+1) ) 
XFREQ(I+1)=FREQ( 1+1 )+l 
9300  CONTINUE 
9200  CONTINUE 
C 

C   CALCULATE  CHI-SQUARE  STATISTIC 
C 

C   ZERO  OUT  XX(I,J),P( I) ,E( I) 
C 

0099  DO  0100  1=1, IN 

pm=o.o 

E( I) =0.0 
DO  0110  J=l,3 
XX( I ,J)=0.0 
0110  CONTINUE 

0100  CONTINUE 
C 

C  LOAD  INTERVAL  BOUNDARI ES ( RANGES )  AND  FREQS  IN  XX  ARRAY 

C  SUBJECT  OT  E(I).GT.5.0. 

C  ALSO  CALCULATE  P(I)  AND  E(I). 

C 

c 

C   CALCULATIONS  FOR  EXPONENTIAL  DISTRIBUTION 
C 

IF(TYDIST.EQ.1)G0  TO  0111 

M=l 

XX(1,1 )=AA 

DO  0200  1  =  1, IN 

TEMP=0.0 

DO  0300  J=M,IN 

TEMP=TEMP+FREQ( J) 

XX(I ,2)=RANGE( J ) 

IF( J.EQ. IN)GO  TO  0320 

XX(I+1»1)=XX(I  2 ) 

P{  I)  =EXP(- (LAMEST* (XX ( 1,1 )-A4) ) )-EXP(-( LAMEST* 
X(XX( I,2)-AA) )) 

E(  I)=N*P(I ) 

IFIE ( I  ).GE.5.0)G0  TO  0310 
0300  CONTINUE 
0310  XX(I,3)=TEMP 

M=J+1 
0200  CONTINUE 
C 

C   FOR  LAST  INTERVAL  TEST  IF  EUJ.LT.5.0,  IF  SO  STORE  IN 
C   PREVIOUS  INTERVAL  AND  RECALCULATE  E(I-l) 
C 
0320  P(I)=EXP(-(LAMEST*(XX( I,1)-AA) ) ) 

E( I)=N*P(I ) 

IF(E(I).GE.5.0)G0  TO  0330 

XX(I-1 ,3)=XX( I-1,3)+TEMP 

XX( I-1,2)=RANGE(  IN) 

XX( I  ,1  )=0.0 

XX( I  ,2)=0.0 

XX( I ,3)=0.0 

P( I-1)=EXP(-(LAMEST*(XX( 1-1,1)- A A)  )  ) 

E( I-1)=N*P(I-1) 

P(  I)  =  0.0 
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E(I)=0.0 
GO  TO  0400 

0330  XX( I ,3)=TEMP 
GO  TO  0400 

C 

C   CALCULATIONS  FOR  HYPEREXPONENT I AL  DISTRIBUTION 
C 
0111  M=l 

XX(1,1)=AA 
DO  0201  1=1, IN 
TEMP=0.0 
DO  0301  J=M,IN 
TEMP=TEMP+FREQ( J ) 
XX(I ,2)=RANGE(J) 
IF( J.EQ.IN)GO  TO  0321 
XX(I+l«L)«XX(If2) 

P( I )=P1*(EXP(-2.0*P1*LAMEST*(XX( I, D-AA) J-EXP (-2 . 0*P 1* 
XLAMEST*(XXU ,2)-AA) ) ) +  P2* ( E XP (-2 .0*P2*LAMEST* 
X(XX( I, D-AA) )-EXP(-2.0*P2*LAMEST*(XX(I f 2)-AA) ) ) 
E( I)=N*P(I) 

IF(E(I ).GE.5.0)GO  TO  0311 
0301  CONTINUE 
0311  XX( I ,3)=TEMP 

M  =  J+1 
0201  CONTINUE 
C 

C   FOR  LAST  INTERVAL  TEST  IF  EUJ.LT.5.0,  IF  SO  STORE  IN 
C   PREVIOUS  INTERVAL  AND  RECALCULATE  E(I-l) 
C 
0321  P( I )=P1*EXP(-2.0*P1*LAMEST*(XX(I , D-AA J ) +P2*EXP (-2.0* 
XP2*LAMEST*(XX(  I,  D-AA)  ) 
Ell) =N*P(I ) 

IF(E( I ).GE.5.0)GO  TO  0331 
XXU-1  ,3)=XX(I-1,3)+TEMP 
XX(I-1 ,2)=RANGE(  IN) 
XX(I, 1)=0.0 
XX(I,2)=0.0 
XX(I ,3)=0.0 

P(  I-U  =  P1*EXP(-2.0*P1*LAMEST*(XX(I-1,D-AA)  )  +  P2* 
XEXP(-2.0*P2*LAMEST*(XX( 1-1,1 )-AA)) 
E(I-1)=N*P(I-U 
P(D=0.0 
E( I)=0.0 
GO  TO  0400 

0331  XX(I,3)=TEMP 
C 

C   CALCULATE  TOTALS  OF  P(I),  E(I)  AND  XX(I,3) 
C 
0400  XSUM=0.0 

PSUM=0.0 

ESUM=0.0 

DO  0500  1=1, IN 

XSUM=XSUM+XX(I ,3) 

PSUM=PSUM+P( I ) 

ESUM=ESUM+E( I ) 
0500  CONTINUE 
C 

C   TEST  TO  FIND  THE  NUMBER  OF  INTERVALS  IN  XX(I,J)  THAT  ARE 
C   IN  ORDER  TO  DETERMINE  THE  NUMBER  OF  DEGRESS  OF  FREEDOM 
C 

DO  0600  1=1, IN 

IF(XX( I ,2) .EQ.O.O)GO  TO  0610 
0600  CONTINUE 
0610  N I  NT  =  1-1 

NDF=NINT-2 
C 
C   CALCULATE  THE  CHI-SQUARE  STATISTIC 


C 


CHISUM=0.0 
DO  0700  I=1,NINT 

CHISUM  =  CHISUM+( ( XX ( I , 3 )-E ( I )  )**2)/E( I) 
0700  CONTINUE 
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c 

C   WRITE  THE  CHI-SQUARE  TEST  RESULTS 
C 

WRITE(6,0502)TITLE 
0502  FORMAT (• 1« ,T10, 'CHI-SQUARE  TEST  RESULTS  F0R:',20A4// 

XT10, 'INTERVALS' ,T2  7, 'FREQ' ,T45, • INTERVAL' ,T65,'EST  OF' 
X/T8,  'LOW  ,T17t  'HIGH' ,T27,' COUNT'  ,T47,'PR0B»  ,T60T 
X'NO.  IN  INTERVAL' ) 
DO  0800  I=1,NINT 

WRITE(6,0602)XX( I, 1) ,XX( 1,2) ,XX( 1,3) ,P( I ) ,E( I ) 
06  02  F0RMAT(F10.2,4X,F6.2,F10.0,F21 -4,F20.2) 
0800  CONTINUE 

WRITE(6,0702)XSUM,PSUM,ESUM,NINT,NDF,CHISUM,XBAR, 
XLAMEST 
0702  FORMAT (T6,  'TOTALS'  ,T21,F10.0,F21.4,F20.2//T10, 

X'NO.  OF  INTERVALS=' , I4//T10, 'DEGREES  OF  FREEDOM= ■  , I  4// 
XT10, 'CHI-SQUARE  ST ATI STI C=« , Fl 0. 2/ /T10 , ■ EST  OF  MEAN=', 
XF10.2//T10,'LAMEST=« ,F10.4) 
IF(INCS.EQ.O)GO  TO  9999 
IN=TEMPIN 
9999  CONTINUE 
STOP 
END 

SUBROUTINE  SORTIA,N) 
DIMENSION  AIN) 
NPASS=Nt-1 
DO  1000  I=1,NPASS 
NSTOP=N-I 
DO  1000  J=1,NST0P 
IFIAl J).LE.A(J+1))G0  TO  1000 
TEMP  =  AU) 
A(J)=A(J+1) 
A( J+1)=TEMP 
1000  CONTINUE 
RETURN 
END 

SUBROUTINE  STAT  I X, N, XBAR, V AR, S ) 
REAL  X(N) 
SUM=0.0 
DO  1  1=1, N 

1  SUM=SUM+X(I) 
XBAR=SUM/N 
VAR=0.0 

DO  2  I=1,N 

2  VAR=VAR+(X(I  )-XBAR)**2 
VAR=VAR/IN-1) 

S=SQRT (VAR) 

RETURN 

END 

SUBROUTINE  HISTOI I N, FREQ ,R ANGE , T ITLE ) 
DIMENSION  FREQ (20) , RANGE (20) , J OUT (20) , TITLE (20) 
DATA  NUTH/'      •/,«/■****•/ 
WRITE  (6,4)  TITLE 

4  FORMAT( • 1* ,///,8X,20A4//) 
JN=IN 

INT  =  0 

IF(JN.GT.O)  GO  TO  10 

INT=1 

JN=-JN 
10  DO  12  1=1, JN 
12  JOUT( I )=FREQ(I ) 

WRITE  (6,5)  (JOUTI I) ,1=1, JN) 

5  FORMAT( 'OFREQUENCY' ,2016) 
WRITE  (6,7) 

7  FCRF,AT(  •  -«  ,128(  '-'  )  ) 
C      FIND  LARGEST  FREQUENCY 
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FMAX=0. 
DO  20  I=1,JN 

IF(FREQU).GT.FMAX)  FMAX=FREQ(IJ 
20  CONTINUE 
SCALE 
JSCAL=1 

IF(FMAX.GT.60.)  JSCAL= ( FMAX+59 . )/60. 
FSCAL=JSCAL 
DO  50  1=1, JN 

50  JOUT( I )=NOTH 
MAX=FMAX/FSCAL 
DO  80  1=1, MAX 
X=MAX-( 1-1) 

DO  70  J=1,JN 
•  IF(FREQ(J)/FSCAL.GE.X)  JOUT(J)=K 
70  CONTINUE 

IX=X*FSCAL 
80  WRITE  (6,2)  IX , ( JOUT ( J ) , J= 1 , JN ) 

2  FORMATt I6,4X,20(2X,A4) ) 
WRITE  (6,7) 

IF( INT.EQ.l)  GO  TO  16 

WRITE  (6,3)  (RANGE(J) ,J=1, JN) 

3  FORMAT (  'OINTERVAL' , 2X, 19 ( F 5. 1 , IX ) , F5 . 1// ) 
GO  TO  15 

16  DO  51  1=1, JN 

51  JOUT( I )=RANGE( I ) 

WRITE  (6,6)  (JOUT( I ) ,I=1,JN) 
6  FORMAT ( 'OINTERVAL'  , 2X , 19 ( I  5 , IX) , 15// ) 
15  WRITE  (6,1)  K,JSCAL 
1  FORMATt 'OEACH  ■ , A4 , '  EQUALS  ',12,'  POINTS'/) 

RETURN 

END 
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TABLES  OF  CHANNEL  SERVICE  STATISTICS  FOR  MAY-JUNF 
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